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Preface 


Volume I 


The second High Speed Research Program Sonic Boom Workshop was 
held at NASA Ames Research Center May 12-14, 1993. The purpose of 
this workshop was to provide a forum for government, industry and 
university participants to present and discuss progress in their research. 
The workshop was organized into three sessions dealing with 
atmospheric propagation, acceptability and configuration design. 
Attendance at the workshop was by invitation only. 

This volume of the workshop proceedings includes papers on 
atmospheric propagation and acceptability studies. Significant progress 
is noted in these areas in the time since the previous workshop a year 
earlier. In particular, several papers demonstrate an improved capability 
to model the effect of atmospheric turbulence on sonic booms. This is a 
key issue in determining the stability and, ultimately, the acceptability, 
of shaped sonic booms. In the area of acceptability, the PLdB metric has 
withstood considerable scrutiny and is validated as a loudness metric for 
a wide variety of sonic boom shapes. The differential loudness of 
asymmetric sonic booms is better understood, too. 


PWeiDfNe PAGE BLANK NOT FH.MID 


v 



N94- 28189 


Progress in Modeling Atmospheric Propagation 

of Sonic Booms 


Allan D. Pierce 

The Pennsylvania State University 
Graduate Program in Acoustics 
157 Hammond Building 
University Park, PA, 16802 


Presentation at NASA HSR Sonic Boom Workshop 
NASA Ames Research Center 
May 12-14, 1993 


i 



The improved simulation of sonic boom propagation through the real atmosphere re- 
quires greater understanding of how the transient acoustic pulses popularly termed 
sonic booms are affected by humidity and turbulence. A realistic atmosphere is in- 
variably somewhat turbulent, and may be characterized by an ambient fluid velocity 
v and sound speed c that vary from point to point. The absolute humidity will also 
vary from point to point, although possibly not as irregularly. What is ideally de- 
sired is a relatively simple scheme for predicting the probable spreads in key sonic 
boom signature parameters. Such parameters could be peak amplitudes, rise times, or 
gross quantities obtainable by signal processing that correlate well with annoyance or 
damage potential. The practical desire for the prediction scheme is that it require a 
relatively small amount of knowledge, possibly of a statistical nature, concerning the 
atmosphere along the propagation path from the aircraft to the ground. The impact 
of such a scheme, if developed, implemented, and verified, would be that it would 
give the persons who make planning decisions a tool for assessing the magnitude of 
environmental problems that might result from any given overflight or sequence of 
overflights. 


Realistic Sonic Boom Propagation Problem 

• Turbulent atmosphere 

• diffraction by smaller turbulent eddies 

• focusing and defocusing by larger turbulent eddies 

• Molecular relaxation important 

• Humidity controls molecular relaxation 

• Nonlinear distortion 

• tendency toward waveform steepening 

• stretching of waveform 

• more rapid elimination of very narrow spikes 

pennSwe 
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The technical approach that has been followed by the author and some of his col- 
leagues is to fonnulate a hierarchy of simple approximate models based on fundamen- 
tal physical principles and then to test these models against existing data. 

For propagation of sonic booms and of other types of acoustic pulses in nonturbu- 
lent model atmospheres, there exists a basic overall theoretical model that has evolved 
as an outgrowth of geometrical acoustics. This theoretical model depicts the sound 
as propagating within ray tubes in a manner analogous to sound in a waveguide of 
slowly varying cross-section. Propagation along the ray lube is quasi-one-dimensional, 
and a wave equation for unidirectional wave propagation is used. A nonlinear term 
is added to this equation to account for nonlinear steepening, and the formulation has 
been carried through to allow for spatially varying sound speed, ambient density, and 
ambient wind velocities. The model intrinsically neglects diffraction, so it cannot take 
into account what has previously been mentioned in the literature as possibly impor- 
tant mechanisms for turbulence-related distortion. The model as originally developed 
could predict an idealized N-wavefonn which often agrees with data in terms of peak 
amplitude and overall positive phase duration. It is possible, moreover, to develop 
simple methods based on the physics of relaxation processes for incorporating molecu- 
lar relaxation into the quasi-one-dimensional model of nonlinear propagation along ray 
tubes. 


Simpler Propagation Models 

• Propagation along ray tubes 

• stratified atmosphere 

• turbulence ignored 

• spiking and roundening effects not predicted 

• Model can account for 

• gross magnification and demagnification 

• nonlinear distortion 

• molecular relaxation contribution to rise times 

pennSkte 
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There are currently two methods that are in use for carrying out such an incorporation 
of relaxation phenomena into propagation predictions; one is a numerical algorithm 
that arose out of the 1973 doctoral dissertation by Pestorius, which carries forward the 
propagation along a ray tube as an alternating sequence of two basic types of steps. 

In one step one has linear propagation of a Fourier superposition of frequency compo- 
nents, and each frequency component is shifted in phase and attenuated in an appro- 
priate manner with propagation along a given distance interval. In the other step, the 
nonlinear distortion is carried out according to inviscid weak shock theory through the 
same distance interval. 

The author and his colleagues, on the other hand, have been working with an explicit 
set of approximate partial differential equations analogous to Burgers' equation, an 
early version of which can be found in the author’s 1981 textbook. One very sim- 
ple model that has been used by the author and his colleagues is what is termed the 
asymptotic quasi-steady theory of sonic boom wavefonns. 


Solution Techniques 


• Transient evolution using Pestorius algorithm 

• split-step algorithm 

• nonlinear distortion step 

• molecular relaxation step using Fourier transforms 

• Asymptotic quasi-steady theory 

• basic waveform shape predicted with neglect of molecular 
relaxation 

• rise-phase corrected for molecular relaxation 

• correction based on local humidity, temperature, and net 
pennState pressure jump in shock 
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The asymptotic quasi-steady theory predicts an explicit waveform shape near the lead- 
ing shock, given the waveform peak amplitude and the local humidity and tempera- 
ture. The model incorporates modecular relaxation, which is slower for dryer air and 
consequently a cause of sharper bangs in humid air. 

This can alternately be termed the “steady-state” model or the “quasi-frozen wave- 
form” model. The terminology is not ideal, and one must first understand the detailed 
assumptions involved before adopting any conceptions about what the terminology im- 
plies. The ideas involved go back to early papers by G. I. Taylor and Richard Becker 
on the structure of shock waves, only here the mechanisms of interest are molecular 
relaxation rather than viscosity or thennal conduction. The first tenet of the theory 
is that molecular relaxation is important only in the rise phase of waveforms. Such 
seems justified because the characteristic times, such as positive phase duration, asso- 
ciated with other portions of the wavefonn are invariably much longer than the char- 
acteristic relaxation times for molecular relaxation. During most of the time at which 
the wavefonn is being received, it is reasonable to assume that the air is in complete 
quasi-static thermodynamic equilibrium. Molecular relaxation is a nonequilibrium ther- 
modynamic phenomenon and is important only when pressure is changing rapidly, 
with characteristic times of the order of a few milliseconds or less. 


Molecular-relaxation correction 

• Rise-phase prediction near t = f . a |,: 

P = (^/dsii/ii.se(^ /. S h, parameters) 

/rise(r, parameters) —» 0 as r — ► -oo 
/rise(T, parameters) -*■ (Ap) sh as r +oo 

• parameters are (Ap) a \ u temperature, and humidity 

• Composite expression: 

V =Pbasic(t)H(t sU - t) 

+ i&phufnsM - < s in parameters) 

+ [pbasic(0 - ( Ap) s u]H{t - t s I,) 

penn State 
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A second hypothesis, which is related to the first, but which requires extensive analy- 
sis for its justification, is that the rise phase of the waveform is detennined solely by 
the peak overpressure of the shock and the local properties of the atmosphere. Strictly 
speaking, one expects the waveform received at a local point to be the result of a 
gradual evolution that took place over the entire propagation path, so it depends in 
principle on the totality of the atmospheric properties along the path. However, an N- 
wave shape, or at least the positive phase portion, is often established fairly close to 
the source (i.e., the flight trajectory in the case of sonic boom generation) relative to 
the overall propagation distance. With increasing propagation distance, the peak over- 
pressure decreases, but does so very slowly, and the positive phase duration increases, 
but also does so very slowly. There is a net loss of energy from the wave and the 
loss takes place almost entirely within the rise phases of the shocks. However, the 
manner in which the peak overpressure decreases and the positive phase duration in- 
creases is virtually independent of the energy loss mechanism. 

One should note in particular that the model based on the asymptotic quasi-steady 
theory predicts rise limes. 



Rise-time prediction 

• (At)risp ' 

to 0.9. 

time interval for / r j se (i-£ s i,, parameters) to rise from 0.1 

• parameters are (Ap) s h, temperature, and humidity 

• (At) r i S e is function of “parameters” 


(AOrise = F([Ap] S h, temperature, humidity) 

where F 

is “universal” function. 

penn State 
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The viewgraph here sketches the principal ideas that are embodied in calculating the 
rise-phase of a sonic boom profile according to the quasi-static theory. The parameter 
£ in the diagram is x — v s | ,t. Ahead of the shock the overpressure p is asymptoti- 
cally zero, and die theory predicts the manner in which it approaches zero. Behind 
the shock p asymptotically approaches die net shock overpressure P s h. and here again 
the theory predicts the maimer in which this asymptotic limit is achieve. For points in 
between one must must numerically solve a set of coupled ordinary nonlinear differen- 
tial equations. One interesting aspect of die solution is that the nitrogen relaxation is 
only important is the later portion of this rise phase. 


Early rise phase: O, relaxation dominates 
Later rise phase: N, relaxation dominates 


The theoretical rise phase is determined using asymptotic 
and numerical solution methods: 
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Kang carried out detailed comparisons of the predictions of the frozen-profile model 
with actual wavefonns of sonic booms, recorded by the US Air Force in the Mojave 
Desert in 1987. The original comparison reported in Kang’s doctoral thesis, unfortu- 
nately, was flawed because the reflection at the ground was incorrectly taken into ac- 
count. (That such may have been the case was first suggested to the author by Gerry 
McAninch as a result of a conversation with Alan Wenzel.) 


Waveforms measured during flight tests 
at ground level 

• (Af) r i. se can be measured for each data sample. 

• Theory predicts 

(Af) r j se = F([Ap] s |,, temperature, humidity) 
where F is “universal” function. 

• For comparison of data with theory, what value of [Ap] s |, does one 
use in the “universal function” F? 

• Theory is based on idealization of plane wave propagating in one 
direction through unbounded medium. Measurements were made 
on the ground 

pennSwe 
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The corrected procedure takes the ground as rigid and the reflection process as lin- 
ear, so that the wavefonn at the ground has the same time dependence as the incident 
wave, only the amplitude is twice as great. The theoretical predictions based on accu- 
mulated nonlinear effects for a unidirectional propagating wave are applicable to the 
so-inferred incident wave. 


Rigid ground idealization: 


Pground — 2p mc , (Ap)sli,grnd — 2(Ap) s h 

»inc 

• Theory predicts a function F, where 

(A t) rise = F([Ap] s |,,i„c, temperature, humidity) 
decreases with increasing (A p) roughly as (Ap) -1 

• Previous comparisons were erroneously based on above with 
argument taken as [Ap] S |, ignK | rather than (l/2)[Ap] s h igriK j 

• Using a (A p) that is too large by a factor of 2 means you tend to 
underpredict the rise time. 

penn State 
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This summarizes the comparison of rise lime data with the asymptotic quasi-steady- 
stale theory. The overpressures on the horizontal axis are those actually observed in 
waveforms recorded at the ground. The theoretical curve is derived assuming that the 
incident wave’s overpressure is one-half of what is measured. All of the data was 
taken at limes when the humidity and temperature were very nearly the same, so one 
theoretical curve suffices for the entire data set. The data was taken in 1987 in the 
Mojave desert, with various airplanes, fiying at various altitudes and with various 
Mach numbers. The relative humidity was 24% and the temperature at the ground 
was 33° C. 



The overall result of the comparison, with ground reflection taken into account as here 
described, is that the theoretically predicted rise times are roughly the same as the av- 
erage rise time of the experimental wavefonns under conditions of the same incident 
waveform pressure amplitude and the same atmospheric humidity. 


Inferences from updated theory-data comparison: 

• Relaxation theory predicts rise times of correct order of magnitude 

• Theoretical predictions of rise times tend (but not in all cases) to 
be lower than observed in field data 

• Turbulence is major factor in rise times. 


pennSiate 
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The rise phase structure of die waveform is basically a tug-of-war between nonlin- 
ear steepening and molecular relaxation. When the boom passes through a region 
where the molecular relaxation is weaker, the nonlinear steepening causes the wave- 
fonn to sharpen up and causes the rise time to decrease until the mechanisms balance 
each other out. One can associate a characteristic adjusUnent time with this restora- 
tion of the balance between the two mechanisms. The quasi-steady hypothesis used 
in the simpler models hypothesis rests on the assertion that this characteristic adjust- 
ment time is substantially less than any characteristic time it takes for the waveform 
to propagate over a path segment within which the relevant atmospheric properties 
(especially the absolute humidity) change appreciably. A current question regarding 
the closely related and competing effects of molecular relaxation and nonlinear steep- 
ening is just how resilient is the steady-stale model. 

Raspet has referred to the characteristic distance over which recovery from a perturba- 
tion to the asymptotic waveform lakes place as the healing distance. 


Concept of healing distance: 


• Suppose rise phase of waveform slightly perturbed from asymptotic 
quasi-steady-state form 

• For further propagation through a homogeneous medium, the 
perturbation dies out 

• Rise phase eventually evolves to asymptotic fonn that depends on 
x and t only in combination x - V^t. 

• What is characteristic additional propagation distance for the 
perturbation to die out? 


pennSwe 


12 




The question of the magnitude of the healing distance has been answered tentatively 
by detailed numerical computations of transient evolution of waveforms over large 
propagation distances by Raspet and others, with the apparent prediction that it takes 
propagation distances of several kilometers (the value depending on the peak ampli- 
tude) for the wavefonn to recover from slight perturbations in the steady-state shape. 
Even more so than is the case for the rise time, there is room for considerable arbi- 
trariness in the definition of this healing distance. The present author suspects that 
one can devise a meaningful definition for which the numerical value of this healing 
distance is less than a kilometer for representative cases of interest. 

That the latter speculation has some credibility can be seen at once when one consid- 
ers that a typical value for the pertinent relaxation time is about 1 ms (correspond- 
ing to the relaxation time of N 2 in air with 50% relative humidity). The waveform 
moves with roughly the sound speed, which is of the order of 340 m/s, and a pressure 
amplitude of 50 Pa would move with an additional speed increment of (5Pj poc = 
(1.2)50/400 = 0.15 m/s. If such a peak lags a zero-crossing by a distance of 
(.001) (340) m, then the distance for it to overtake the zero-crossing in the absence 
of any dissipation effect is approximately (.001)(340) 2 /0.15 or 0.8 km. To put such 
an estimate in perspective, one can contrast this with a distance of 1 1 km for a typical 
height of the tropopause and with a representative distance of 15 km for a ray trajec- 
tory from the aircraft flight track to the ground. The numbers sometimes mentioned 
for the thickness of the atmosphere’s turbulent boundary layer, on the other hand, are 
much smaller, on the order of 1 to 2 km. 


Order of magnitude of healing distance: 


• Take pertinent relaxation time as 1 ms 

• Waveform moves with speed c. « 340 m/s 

• Pressure amplitude of 50 Pa moves with additional speed increment 
of 

fi{Ap)/p 0 c = (1.2)(50)/(400) = 0.15m/s 

• Propagation distance for peak to overtake zero crossing in absence 
of dissipation: 

(0.001 )(340) 2 /0. 15 « 0.8km 

PENNStATE 
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Kang and the author some time ago initiated a systematic study of the tendency to- 
ward the steady-state profile in which the analysis was based more on analytic con- 
siderations rather than lengthy numerical case studies. In this theory the steady-state 
profile provides a framework whereby perturbations to the profile can be regarded as 
a superposition of natural eigenfunctions of fixed shape as seen by someone moving 
with the nominal shock speed. Each such natural eigenfunction has its own natural 
decay time which results as an eigenvalue in the theory. The task then emerges of 
systematically determining these eigenfunctions and the associated eigenvalues. 


Eigenfunctions and characteristic healing distances: 


• Propagation equations are 

du du du / d^u v — > du„ 

m +c dx + liu !h'. ~ 6c 'd + ^ !h = 0 

du u 1 du 

IT ^ U{J = 

at Tj/ at 

where p ~ poxu . 

• Steady state solution is 

u = F((); u„ = F„(0 

where £ = x - Ksh t. 

• Take perturbed solution to be of form 

pennSwe 
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This summarizes the mathematical problem that ideally should be solved to determine 
a sequence of healing distances. One expects the so-posed eigenvalue problem to have 
several solutions that correspond to real eigenfunctions. The structure of the problem 
is still to be studied, and one does not have any orthogonality theorems as yet regard- 
ing different eigenvector functions. However, a crude solution can be found if one 
replaces the governing equations by Burgers’ equation with an effective bulk viscosity. 


Equations governing healing eigenfunctions: 

tt = F{£) + iP(£)e~ Xx -, Ur = Frit) + MOe~ Xx 

are inserted into propagation equations; one keeps only linear 
terms in the ip's, with result in matrix form 

[£]M = A [M]{ip} + X 2 [M]{ip} 

where [£] , [ M ] , and [TV] are 3-by-3 matrices made up of linear 
operators, each possibly involving differentiation with respect to £ 
and the steady state profile functions F(£), F^{£), F2{£) 

• {ip} is an eigenfunction array (ip,ip\,ip 2 ) 

• Nontrivial solution (boundary condition of ip ' s equal to zero at 
£ = ±oo) exists only for special values of A. 

• Special values of A construed as eigenvalues and as reciprocals of 
p^^^haracteristic healing distances. 
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Predictions of healing distance 
based on effective bulk viscosity: 


Tisza’s equivalence: 


/ib»ik » 2pe(Ac) I/ r I , 


• so we infer 


^lieal = 16 cr„ 


.2 \2 


pc 

/j(A p) sh 


pennState 


Representative numerical values 
based on effective bulk viscosity: 


f'lieal = lGCT,, 


c ) sli 


• Typical numerical values (0 2 relaxation) 


(Ap) 8 h = 50 Pa; . • = 2300 

/j(Ap) s |, 


= 3 x 10 -4 ; Tjy = 10 -5 s 


• so we infer L|, e ai ~ 250 in 

• But you can get much larger values using parameters for N 2 
relaxation. 


pennState 
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The author and his colleagues have recently been exploring various methods for com- 
putation of sonic boom propagation through turbulent aUnospheres and have obtained 
a generalization of the Burgers equation which has some similarities to the KZK equa- 
tion and to the NPE equation of McDonald and Kuperman. The equation can be re- 
garded as a string of “small terms” lacked onto the inviscid linear Burgers equation, 
with individual tenns accounting for nonlinear steepening, viscous attenuation, refrac- 
tion, molecular relaxation, and diffraction. 


The Penn State Univ Propagation Equation (PSUPE): 

• Generalization of Burgers equation (which really should be called 
the Cole equation, as inferred from Cal Tech literature of late 
1940’s by ADP) 

• Term for diffraction by smaller turbulent eddies 

• Molecular relaxation term 

• Nonlinear steepening term 

• Turbulence can be simulated using Fourier transforms (series) 
with random number algorithms used in selection of coefficients. 

• Larger scale turbulence and ambient atmospheric stratification can 
be incorporated in multiplicative Blokhintzev factor (roughly same 
as 1/ yj 'A, where A is ray tube area) which varies with distance 
along central ray. 

penn State 
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Pestorius had a good idea 


• Basic Pestorius algorithm alternated between nonlinear distortion 
(NL) and absorption-dispersion (AD) steps 

• Noise-Con 93 paper by ADP shows that this is rigorously correct 
in limit of small step sizes Ax. Taking limit yields partial differ- 
ential equation 

du 1 du .. , , . . , , 

& + “-^nlM+AWu} 

where right-side terms are same as appear in PSUPE. 

• Why not use same idea to handle the diffraction term? 
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Concluding remarks 

• Solving PSUPE to realistically simulate sonic boom statistics will 
require major theoretical innovations in computational acoustics. 
(But so what?) 

• Two doctoral theses presently in progress at Penn State on 
alternate approaches to solving PSUPE 

• Kirchhoff integral for the diffraction step. 

• Finite-difference algorithms using flux-corrected transport 
methods of Boris and McDonald. 

• Another idea being pursued by the present author indepen- 
dently is that the healing eigenfunctions are natural basis set 
for a variational (or Galerkin) formulation which reduces the 
dimensionality of the problem, and which leads to simpler ways 
of decomposing messy waveforms encountered in field data. 


pennSwe 
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Good morning. My name is Dr. Victor W. Sparrow of the Penn State 
University Graduate Program in Acoustics. The Graduate Program in 
Acoustics is a Department in the Penn State College of Engineering. My 
co-author is Tom Gionfriddo, a graduate student at Penn State who 
finished up his M.S. degree in Acoustics early in the Fall of 1992. Much 
of the work I will be presenting today is the result of Tom's effort on his 
master's thesis. 

The topic I will be discussing today is Implications for High Speed 
Research: The Relationship Between Sonic Boom Signature Distortion 
and Atmospheric Turbulence. But before we get to these implications, 
let us review a little history concerning previous research on sonic 
boom waveform distortion. 


Implications for High Speed Research: 

The Relationship Between Sonic Boom Signature Distortion 
and Atmospheric Turbulence 


Victor W. Sparrow and Thomas A. Gionfriddo 
Graduate Program in Acoustics 
The Pennsylvania State University 


Presentation at NASA HSR Sonic Boom Workshop 
NASA Ames Research Center 
May 12-14, 1993 
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In 1968 Dr. Allan Pierce hypothesized that the cause of sonic boom 
distortion, which takes the form of spiked or rounded waveforms, was 
due to atmospheric turbulence. This was a theoretical result, and was 
not widely accepted at the time due to the lack of experimental evidence. 

In 1973 Ribner, Morris, and Chu performed laboratory experiments 
which showed that one could cause sonic boom shaped waves to spike or 
become rounded, if the waves were propagated through a turbulent jet. 
This laboratory result gave some evidence that turbulence could, in fact, 
be the cause of sonic boom waveform distortion in the atmosphere. 
Others also performed similar laboratory experiments. 

In the mid to late 1970's, however, the role of molecular relaxation 
absorption in sonic boom propagation had not yet been established. The 
relative importance of molecular relaxation and atmospheric 
turbulence for sonic boom distortion was not clear. 


• In 1968 Pierce hypothesized (Ref. 1) that the cause of sonic boom 
distortions, such as spiked or rounded waveforms, was due to 
atmospheric turbulence. 

• In 1973 Ribner, Morris, and Chu found in the laboratory (Ref. 2) 
that sonic boom shaped acoustic waves indeed were distorted in a 
turbulent jet, producing both spiked and rounded waveforms. 

• However, the relative importance of atmospheric turbulence and 
molecular relaxation effects had yet to be established. 
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By the early 1980's the theory for molecular relaxation absorption In 
the atmosphere was fairly well understood. The two dominant process 
are Oxygen and Nitrogen relaxation, with humidity (water vapor) being 
the next most important process. The result of the theory is that 
molecular relaxation cannot cause the spikes on sonic boom waves, 
although they can round sonic boom waves somewhat. It is thought that 
the rounding effect is insufficient to explain observed distorted 
waveforms, however. 

After the molecular relaxation theory was understood, it became the 
common notion that atmospheric turbulence is primarily responsible 
for sonic boom distortion. This is an assumption which most workers in 
sonic boom propagation have adopted, since the spiking and rounding 
could not be due to molecular relaxation. Most of the talks during the 
rest of this session make this assumption. 


• By the early 1980’s molecular relaxation was fairly well understood. 
Molecular relaxation cannot cause spikes on sonic boom waves. 

• Thus, most researchers in sonic boom propagation have assumed 
that turbulence must be responsible for sonic boom distortion. 

• Most of the talks in this session make such an assumption. 


pennStot 
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Until recently this assumption has not been tested statistically. Such a 
test would provide a firm foundation for much of the ongoing work on 
sonic boom propagation through turbulence at a number of NASA 
Contractor sites, including The University of Mississippi, The University 
of Texas at Austin, Penn State University, Wyle Laboratories, etc. 

One supposes here that an originally undistorted sonic boom 
propagating through turbulence should, on average, be more distorted 
as it propagates through more turbulence. 


However, 

until recently this assumption has not been tested statistically with 
real sonic boom data and real atmospheric turbulence. 

One supposes that a sonic boom propagating through more turbulence 
should, on average, be more distorted. 


pennSkte 
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Therefore, the purpose our research study is was to test the above 
hypothesis rigorously. That is, the specific purpose is to see if 
increasing travel distances through turbulence is correlated with 
increasing sonic boom wave distortion. This paper documents the 
results of our study. 

In this study it is assumed that the strength of the atmospheric 
turbulence is somewhat uniform, and it is the travel distance of booms 
through the turbulence that is important. This assumption is necessary 
due to the absence of direct turbulence measurements to complement 
the sonic boom experimental data which will be used to test the 
hypothesis. 


The Purpose of this study is to test the above hypothesis rigorously. 


More specifically, is it true that 


the further a boom travels through turbulence 


increased waveform distortion 
? 


pennSkie 
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In this study there were two primary tasks. The first was to develop an 
algorithm for quantifying the distortion in a sonic boom. Such an 
algorithm should be somewhat automatic, with minimal human 
intervention. Once the algorithm was developed, it was used to test the 
previously mentioned hypothesis. This hypothesis testing was the 
second task. Using readily available sonic boom data, we statistically 
tested whether there was a correlation between the sonic boom 
distortion and the distance a boom traveled through atmospheric 
turbulence. 


In this study we 

A. Developed an algorithm to quantify the distortion in a sonic boom 
waveform. 

B. Tested the correlation between this distortion and the distance a 
boom traveled through atmospheric turbulence. 
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The terminology that is used in our paper is described here. The booms 
have a maximum shock overpressure after some rise time. This 
maximum shock overpressure is called the bow shock. The duration is 
then defined as the time as waveform slopes off to the minimum shock 
overpressure at the tail shock. For most of the booms examined, the 
duration was between 75 and 200 milliseconds, and the amplitudes varied 
between 30 and 200 pascals. Most booms have a subsonic fundamental 
frequency in the range of 6 to 10 hertz. Our definition of rise time is 
from 10% to 90% of the maximum shock overpressure. 


SONIC BOOM WAVEFORM TERMINOLOGY 



Duration typically 75-200 ms, amplitude 30-200 Pa 

Subsonic fundamental - 6-10 Hz 

Rise time is time from 10% to 90% of maximum 
shock overpressure 

Why annoying? Rise phase structure important 
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Sonic booms can be distorted in many ways as they propagate through 
the atmosphere. Here a large number of sonic booms were collected 
into categories, and a representative waveform example is shown from 
each category. The waves which were the most undistorted were called 
Classic N. Waveforms showing one large peak were called Peaked. 

Many of waveforms have two distinct peaks, and were called Double- 
peaked. Some waveforms had many peaks, and these were called Multi- 
peaked. The U-wave category was defined as those waves having very 
large spikes on both the bow and tall shocks, the spikes dominating all 
features. The Rounded waveform category had rounded bow and tail 
shocks. All other waveforms, which could not be classified in one of the 
previous categories, were called Messy (for lack of a better term). 


DATA CLASSIFICATION BY 
WAVEFORM SHAPE 



Classic N 



Peaked 





Hh tim. — ■i L— 

> 6su 

Rounded 
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Because of the great variability in the distortions a sonic boom wave can 
undergo, an automatic numerical procedure was developed for 
quantifying the distortion in a sonic boom wave. The first component 
of this quantification is to find a basis for comparison. Since one 
usually puts microphones on the ground, and measures the sonic boom 
wave only after it has been distorted, it is necessary to estimate the wave 
shape of the sonic boom before it was distorted. 

The assumption made here is that before any waveform distortions 
occurred that the sonic boom wave had the shape of a perfect N-wave 
with zero rise times on both bow and tail shocks. The energy in the 
distorted sonic boom is measured, and then it is assumed that the 
undistorted ideal N-wave has the same energy. Obviously, this is an 
approximation. 

Given this information the proper maximum overpressure, duration, 
and start time offset of the ideal N-wave is automatically computed. 
Additional details on the elaborate algorithm used to calculate the 
parameters for the ideal N-wave, given the distorted sonic boom wave, 
are available in the M.S. thesis of Gionfriddo (Ref. 3). 


WAVEFORM DISTORTION ANALYSIS: 
COMPARE DATA TO IDEAL N-WAVE MODEL 



Model of the sonic boom signature just prior to entering the TBL. 

Ideal N-wave and the recorded data have equal acoustic energy. 

The ideal N-wave is superimposed on the same time axis as the recorded data. 


The proper maximum overpressure, duration, and start time offset of the 
ideal N-wave must be determined. 

It is desired to have Classic N data with the lowest mean- squared deviation for any 
waveshape. Therefore, Classic N data are used as a reference for superimposing the 
ideal N-wave correctly over the recorded data. 
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Here is a typical recorded sonic boom wave with its ideal N-wave 
superimposed. 

To quantify the distortion in the measured waveform, the notion of a 
mean-squared deviation is used, defined below. The functions pidealtn] 
and Precordedtn] are both assumed to be digitized data. For the example 
waveform shown here, the mean-squared deviation is 0.075. 


EXAMPLE OF RECORDED BOOM WITH 
IDEAL N-WAVE SUPERIMPOSED 



MEASURE OF WAVEFORM DISTORTION: 
MEAN-SQUARED DEVIATION 


The mean-squared deviation (MSD) 
is defined as: 


MSD s 


2 ( Pidealfn] - Precordedtn] ) 
X (Pidealfn]) 2 
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Here are two more sonic boom waveforms with superimposed ideal N- 
waves. The upper waveform is a Peaked sonic boom wave, and it shows a 
mean-squared deviation of 0.095. The lower waveform is a Rounded 
waveform, having a much larger mean-squared deviation of 0.23. 




-50 0 50 100 150 200 

Time (ms) 


30 




The particular sonic boom data we analyzed was taken near Edwards Air 
Force Base in the late summer of 1987 by the U. S. Air Force. 
Autonomous Boom Event Analyzer Recorder (BEAR) systems took data 
over several days from a wide variety of supersonic aircraft: F-4, F-14, 
F-15, F-16, F-18, SR-71, T-38, AT-38, and F-l 1 ID. The recorders were 
placed at the mile markers along a road in the area. The aircraft were to 
fly perpendicular to the road over a specific flight track. From 44 
aircraft flights, over 500 data files were obtained for analysis. The 
specific position of where the aircraft overflew the road was recorded, 
and this information has been taken into account in our analysis. 


MOJAVE DESERT SONIC BOOM ACQUISITION SCHEMATIC 

F-4, F-14, F-15, F-16, F-18, SR-71, T-38, 

AT-38, and F-l 1 ID aircraft 

Lateral array of thirteen Boom Event 
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To determine the path length a sonic boom will traverse through the 
turbulent boundary layer near the earth's surface, a short exercise in 
three-dimensional solid geometry is needed. Knowing the altitude of 
the aircraft, its Mach number, the lateral ground distance of the 
receiving microphone from the aircraft’s actual flight track, and the 
thickness of the boundary layer, this path length can be obtained. The 
path length is shown as a dark solid line in the diagram. 



Function of: 
altitude 
Mach number 
lateral ground distance 
TBL thickness 
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The turbulent boundary layer is defined as the thickness of the mixing 
layer of the planetary boundary layer. One may assume that the 
turbulence In the mixing layer is somewhat evenly distributed and 
homogeneous. 

To determine the thickness of this layer, a numerical model by A. K. 
Blackadar was employed. The numerical model takes into account 
information from rawlnsonde launches, surface weather data, satellite 
cloud photos, and soil parameters for the site of the sonic boom tests. 
Blackadar's model provides daily profiles to 2000 m height for 
temperature, water content, wind, and boundary layer thickness. 


PATH LENGTH THROUGH THE 
TURBULENT BOUNDARY 
LAYER (CONTINUED) 



TBL thickness estimated using numerical model by A. K. Blackadar 
(PSU meteorology) 

Input atmospheric information from rawinsonde launches, surface 
weather station, satellite cloud photos, soil parameters 

Model provides diurnal profiles to 2000 m for temperature, water 
content, wind, and TBL thickness 
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For the days of the tests, 3 - 8 August 1987, profiles of the turbulent 
boundary layer thickness were obtained from the Blackadar model. One 
can see the thickness of the boundary layer generally grew during the 
day between 7:00 AM and 2:00 PM. Because of the meteorological 
conditions present, the boundary layer grew much more on 3 August 
and 4 August than it did on the other days of the tests. 


TURBULENT BOUNDARY LAYER THICKNESS 
ESTIMATION USING BLACKADAR MODEL 


Local time 


05.00 AM 0700 AM 0900 AM 1100 AM 0100 PM 0300 PM 



1 

M 

*5 
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Before continuing, the analysis procedure will be summarized. On one 
track the BEAR sonic boom data was obtained and calibrated. From this 
data pressure versus time plots were obtained, which were subsequently 
sorted into waveshape categories. Given these plots, the sonic boom 
distortion quantification algorithm was run, and mean-squared 
deviations from the computer generated corresponding ideal N-waves 
were obtained for all the the waveforms. 

On the other track, the aircraft flight parameters and geometry were 
combined with the meteorological data and subsequent predictions of 
the boundary layer height from Blackadar's model. From this 
information the path length through the turbulence was found for 
each recorded sonic boom waveform. 

It was then possible to determine if a statistical correlation existed 
between the mean-squared deviation and the path length the sonic 
boom traveled through the turbulence. 


ANALYSIS PROCEDURE 
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During the analysis procedure it became Immediately apparent that 
those waveforms which were either shaped as a Classic N-waves, or 
were nearly shaped as such, were primarily manifest only in the early 
morning hours. 

This result leads us to believe that as the turbulent boundary layer grew 
through the day, that the number of undistorted waveforms decreased. 
This result is averaging over all of the usable observation data. 


SIMPLE DEMONSTRATION: FRACTION OF RECORDED 
SONIC BOOM WAVEFORMS THAT ARE LEAST-DISTORTED 



8. ® E E ■ i. 3.^— 

333312 


t- so 9 o 


S3 


445 recordings from 7/31-8/7 
overflights during 7AM- 1PM 

Fraction is number of classic N 
and nearly classic N waveforms 
per number of waveforms 
recorded during the hour intervals 

classic N ^ least distorted 


number of 

recorded 

waveforms 


TBL thkknees Increasing 


36 


This is a plot of the mean-squared deviation as a function of altitude (or 
roughly Mach number, since faster planes generally flew higher). It is 
apparent that increased altitude and speed imply decreased waveform 
distortion and a decreased spread of data points. Now higher and faster 
flying planes generally will have shorter propagation paths through 
the turbulent boundary layer, which can be shown from simple 
geometry. Thus, it appears as if longer propagation paths through the 
turbulence result in larger mean-squared deviations, i.e., more distorted 
waveforms. 


ALTITUDE AND MACH NUMBER INFLUENCE UPON 
WAVEFORM DISTORTION 


Increased altitude and speed decreased waveform distortion and 

decreased spread of data points 
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Here are some results listed by plane type of the linear correlation 
coefficients between the mean-squared deviation and the path length 
through the turbulent boundary layer. It is seen that there is strong 
correlation in some cases ( F-18 and F-15) and fair correlation in the 
others. This also gives us some evidence that increasing distortion is 
correlated with increasing path length through the turbulence. 

The F-4 and F-16 data have the lowest correlations and the lowest 
altitudes. For these cases the sonic boom signatures probably did not 
have time to develop into an N-wave by the time it began to interact 
with the turbulence. The other aircraft types flown had too few data 
points to draw any statistical conclusions. 

We also are currently working on obtaining correlation coefficients 
grouped by altitude and mach number as well as by plane type, to 
determine how these factors Interrelate. 


CORRELATION COEFFICIENTS FOR 
MEAN-SQUARED DEVIATION AND 
TBL PATH LENGTH 


Aircraft type 

Correlation coefficient: MSD vs 
TBL path length 

Number of samples 

F-18 

0.712 

62 

F-15 

0.591 

75 

SR-71 

0.398 

48 

F-4 

0.324 

46 

F-16 

0.318 

65 


(Linear regression correlation coefficients) 


Strong correlation in some cases, fair correlation in others 

F-4 and F-16 have lowest correlations and the lowest altitudes 

- Sonic boom signature before TBL probably not N-wave 

- Less distance for nonlinear steepening to work before TBL 
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The conclusions of this study are the following: a strong linear 
correlation exists between the mean-squared deviation and path length 
through the turbulence for the F-18 and F-15 sonic boom data. Fair 
correlation exists for the SR-71, F-4, and F-16 data. An increase in 
altitude and speed results in decreased waveform distortion and a small 
deviation between distortion values. Looking at the waveform 
classification results, the large percentage of Classic N-wave data 
during each day's early flights seems to correspond with the thin 
boundary layer at that time. 


Conclusions: 

• A strong linear correlation exists between mean-squared deviation 
and path length through the turbulence for F-18 and F-15 sonic 
boom data. 

• Fair correlation exists for SR-71, F-4, and F-16 data. 

• An increase in altitude and speed results in decreased waveform 
distortion and a smaller deviation between distortion values. 

• Looking at waveform classification results, the large percentage of 
classic N wave data during early flights seems to correspond with 

the thin boundary layer at that time. 
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The implications of this study for high speed research are the 
following: Increased interaction between real atmospheric turbulence 
and actual sonic boom data does imply more distorted waveforms. The 
common assumption prevailing in the sonic boom propagation 
community for the last several years has been validated statistically. 
And most importantly, it is now clear that atmospheric turbulence will 
determine how well a shaped sonic boom will remain shaped as it 
propagates to the ground. We are now led to believe that higher and 
faster aircraft having shaped sonic booms will, on average, have more 
shaped boom preserved than will aircraft flying at lower altitudes and 
slower speeds, since flying higher and faster minimizes the path length 
through the turbulence. 


Implications for High Speed Research: 

• Increased interaction between real atmospheric turbulence and 
actual sonic booms does imply more distorted waveforms. 

• This common assumption has been validated statistically. 

• Most Importantly: Atmospheric turbulence primarily will 
determine how well a shaped sonic boom will remain shaped as 
it propagates to the ground. 
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Good morning. My name is Dr. Victor W. Sparrow of the Penn State 
University Graduate Program in Acoustics. The Graduate Program in 
Acoustics is a Department in the Penn State College of Engineering. My 
co-author is Tom Gionfriddo, a graduate student at Penn State who 
finished up his M.S. degree in Acoustics early in the Fall of 1992. Much 
of the work I will be presenting today is the result of Tom's effort on his 
master's thesis. 

The topic I will be discussing today is Implications for High Speed 
Research: The Relationship Between Sonic Boom Signature Distortion 
and Atmospheric Turbulence. But before we get to these implications, 
let us review a little history concerning previous research on sonic 
boom waveform distortion. 


Implications for High Speed Research: 

The Relationship Between Sonic Boom Signature Distortion 
and Atmospheric Turbulence 


Victor W. Sparrow and Thomas A. Gionfriddo 
Graduate Program in Acoustics 
The Pennsylvania State University 


Presentation at NASA HSR Sonic Boom Workshop 
NASA Ames Research Center 
May 12-14, 1993 
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In 1968 Dr. Allan Pierce hypothesized that the cause of sonic boom 
distortion, which takes the form of spiked or rounded waveforms, was 
due to atmospheric turbulence. This was a theoretical result, and was 
not widely accepted at the time due to the lack of experimental evidence. 

In 1973 Ribner, Morris, and Chu performed laboratory experiments 
which showed that one could cause sonic boom shaped waves to spike or 
become rounded, if the waves were propagated through a turbulent jet. 
This laboratory result gave some evidence that turbulence could, in fact, 
be the cause of sonic boom waveform distortion in the atmosphere. 
Others also performed similar laboratory experiments. 

In the mid to late 1970’s, however, the role of molecular relaxation 
absorption in sonic boom propagation had not yet been established. The 
relative Importance of molecular relaxation and atmospheric 
turbulence for sonic boom distortion was not clear. 


• In 1968 Pierce hypothesized (Ref. 1) that the cause of sonic boom 
distortions, such as spiked or rounded waveforms, was due to 
atmospheric turbulence. 

• In 1973 Ribner, Morris, and Chu found in the laboratory (Ref. 2) 
that sonic boom shaped acoustic waves indeed were distorted in a 
turbulent jet, producing both spiked and rounded waveforms. 

• However, the relative importance of atmospheric turbulence and 
molecular relaxation effects had yet to be established. 
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By the early 1980's the theory for molecular relaxation absorption in 
the atmosphere was fairly well understood. The two dominant process 
are Oxygen and Nitrogen relaxation, with humidity (water vapor) being 
the next most important process. The result of the theory is that 
molecular relaxation cannot cause the spikes on sonic boom waves, 
although they can round sonic boom waves somewhat. It is thought that 
the rounding effect is insufficient to explain observed distorted 
waveforms, however. 

After the molecular relaxation theory was understood, it became the 
common notion that atmospheric turbulence is primarily responsible 
for sonic boom distortion. This is an assumption which most workers in 
sonic boom propagation have adopted, since the spiking and rounding 
could not be due to molecular relaxation. Most of the talks during the 
rest of this session make this assumption. 


• By the early 1980’s molecular relaxation was fairly well understood. 
Molecular relaxation cannot cause spikes on sonic boom waves. 

• Thus, most researchers in sonic boom propagation have assumed 
that turbulence must be responsible for sonic boom distortion. 

• Most of the talks in this session make such an assumption. 
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Until recently this assumption has not been tested statistically. Such a 
test would provide a firm foundation for much of the ongoing work on 
sonic boom propagation through turbulence at a number of NASA 
Contractor sites, including The University of Mississippi, The University 
of Texas at Austin, Penn State University, Wyle Laboratories, etc. 

One supposes here that an originally undistorted sonic boom 
propagating through turbulence should, on average, be more distorted 
as it propagates through more turbulence. 


However, 

until recently this assumption has not been tested statistically with 
real sonic boom data and real atmospheric turbulence. 

One supposes that a sonic boom propagating through more turbulence 
should, on average, be more distorted. 
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Therefore, the purpose our research study is was to test the above 
hypothesis rigorously. That is, the specific purpose is to see if 
increasing travel distances through turbulence is correlated with 
increasing sonic boom wave distortion. This paper documents the 
results of our study. 

In this study it is assumed that the strength of the atmospheric 
turbulence is somewhat uniform, and it is the travel distance of booms 
through the turbulence that is important. This assumption is necessary 
due to the absence of direct turbulence measurements to complement 
the sonic boom experimental data which will be used to test the 
hypothesis. 


The Purpose of this study is to test the above hypothesis rigorously. 


More specifically, is it true that 


the further a boom travels through turbulence 


increased waveform distortion 
? 
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In this study there were two primary tasks. The first was to develop an 
algorithm for quantifying the distortion in a sonic boom. Such an 
algorithm should be somewhat automatic, with minimal human 
intervention. Once the algorithm was developed, it was used to test the 
previously mentioned hypothesis. This hypothesis testing was the 
second task. Using readily available sonic boom data, we statistically 
tested whether there was a correlation between the sonic boom 
distortion and the distance a boom traveled through atmospheric 
turbulence. 


In this study we 

A. Developed an algorithm to quantify the distortion in a sonic boom 
waveform. 

B. Tested the correlation between this distortion and the distance a 
boom traveled through atmospheric turbulence. 
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The terminology that is used in our paper is described here. The booms 
have a maximum shock overpressure after some rise time. This 
maximum shock overpressure is called the bow shock. The duration is 
then defined as the time as waveform slopes off to the minimum shock 
overpressure at the tail shock. For most of the booms examined, the 
duration was between 75 and 200 milliseconds, and the amplitudes varied 
between 30 and 200 pascals. Most booms have a subsonic fundamental 
frequency in the range of 6 to 10 hertz. Our definition of rise time is 
from 10% to 90% of the maximum shock overpressure. 


SONIC BOOM WAVEFORM TERMINOLOGY 



Duration typically 75-200 ms, amplitude 30-200 Pa 

Subsonic fundamental - 6-10 Hz 

Rise time is time from 10% to 90% of maximum 
shock overpressure 

Why annoying? Rise phase structure important 
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Sonic booms can be distorted in many ways as they propagate through 
the atmosphere. Here a large number of sonic booms were collected 
Into categories, and a representative waveform example is shown from 
each category. The waves which were the most undistorted were called 
Classic N. Waveforms showing one large peak were called Peaked. 

Many of waveforms have two distinct peaks, and were called Double- 
peaked. Some waveforms had many peaks, and these were called Multi- 
peaked. The U-wave category was defined as those waves having very 
large spikes on both the bow and tall shocks, the spikes dominating all 
features. The Rounded waveform category had rounded bow and tail 
shocks. All other waveforms, which could not be classified in one of the 
previous categories, were called Messy (for lack of a better term). 


DATA CLASSIFICATION BY 
WAVEFORM SHAPE 





U-wave 
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Because of the great variability in the distortions a sonic boom wave can 
undergo, an automatic numerical procedure was developed for 
quantifying the distortion in a sonic boom wave. The first component 
of this quantification is to find a basis for comparison. Since one 
usually puts microphones on the ground, and measures the sonic boom 
wave only after it has been distorted, it is necessary to estimate the wave 
shape of the sonic boom before it was distorted. 

The assumption made here is that before any waveform distortions 
occurred that the sonic boom wave had the shape of a perfect N-wave 
with zero rise times on both bow and tail shocks. The energy in the 
distorted sonic boom is measured, and then it is assumed that the 
undistorted ideal N-wave has the same energy. Obviously, this is an 
approximation. 

Given this information the proper maximum overpressure, duration, 
and start time offset of the ideal N-wave is automatically computed. 
Additional details on the elaborate algorithm used to calculate the 
parameters for the ideal N-wave, given the distorted sonic boom wave, 
are available in the M.S. thesis of Gionfriddo (Ref. 3). 


WAVEFORM DISTORTION ANALYSIS: 
COMPARE DATA TO IDEAL N-WAVE MODEL 



Model of the sonic boom signature just prior to entering the TBL. 

Ideal N-wave and the recorded data have equal acoustic energy. 

The ideal N-wave is superimposed on the same time axis as the recorded data. 

The proper maximum overpressure, duration, and start time offset of the 
ideal N-wave must be determined. 

It is desired to have Classic N data with the lowest mean- squared deviation for any 
waveshape. Therefore, Classic N data are used as a reference for superimposing the 
ideal N-wave correctly over the recorded data. 
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Here is a typical recorded sonic boom wave with its ideal N-wave 
superimposed. 

To quantify the distortion in the measured waveform, the notion of a 
mean-squared deviation is used, defined below. The functions pidealM 
and Precorded[n] are both assumed to be digitized data. For the example 
waveform shown here, the mean-squared deviation is 0.075. 


EXAMPLE OF RECORDED BOOM WITH 
IDEAL N-WAVE SUPERIMPOSED 



MEASURE OF WAVEFORM DISTORTION: 
MEAN-SQUARED DEVIATION 


The mean-squared deviation (MSD) 
ir defined as: 


MSD : 


y, ( Puiealll' J - Precorded [n]) 
X (pideal[nj ) 2 
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Here are two more sonic boom waveforms with superimposed ideal N- 
waves. The upper waveform is a Peaked sonic boom wave, and it shows a 
mean-squared deviation of 0.095. The lower waveform is a Rounded 
waveform, having a much larger mean-squared deviation of 0.23. 



SUPERIMPOSED: 
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The particular sonic boom data we analyzed was taken near Edwards Air 
Force Base in the late summer of 1987 by the U. S. Air Force. 
Autonomous Boom Event Analyzer Recorder (BEAR) systems took data 
over several days from a wide variety of supersonic aircraft: F-4, F-14, 
F-15, F-16, F-18, SR-71, T-38, AT-38, and F-111D. The recorders were 
placed at the mile markers along a road in the area. The aircraft were to 
fly perpendicular to the road over a specific flight track. From 44 
aircraft flights, over 500 data files were obtained for analysis. The 
specific position of where the aircraft overflew the road was recorded, 
and this information has been taken into account in our analysis. 


MOJAVE DESERT SONIC BOOM ACQUISITION SCHEMATIC 

F-4, F-14, F-15, F-16, F-18, SR-71, T-38, 

AT-38, and F-ll ID aircraft 


Lateral array of thirteen Boom Event 
Analyzer Recorder (BEAR) systems 
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To determine the path length a sonic boom will traverse through the 
turbulent boundary layer near the earth's surface, a short exercise in 
three-dimensional solid geometry is needed. Knowing the altitude of 
the aircraft, its Mach number, the lateral ground distance of the 
receiving microphone from the aircraft's actual flight track, and the 
thickness of the boundary layer, this path length can be obtained. The 
path length is shown as a dark solid line in the diagram. 



Function of: 
altitude 
Mach number 
lateral ground distance 
TBL thickness 
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The turbulent boundary layer is defined as the thickness of the mixing 
layer of the planetary boundary layer. One may assume that the 
turbulence in the mixing layer is somewhat evenly distributed and 
homogeneous. 

To determine the thickness of this layer, a numerical model by A. K. 
Blackadar was employed. The numerical model takes into account 
information from rawinsonde launches, surface weather data, satellite 
cloud photos, and soil parameters for the site of the sonic boom tests. 
Blackadar's model provides daily profiles to 2000 m height for 
temperature, water content, wind, and boundary layer thickness. 


PATH LENGTH THROUGH THE 
TURBULENT BOUNDARY 
LAYER (CONTINUED) 



turbulent 
boundary layer 


TBL thickness estimated using numerical model by A. K. Blackadar 
(PSU meteorology) 

Input atmospheric information from rawinsonde launches, surface 
weather station, satellite cloud photos, soil parameters 

Mc^ ™ vides diurnal profiles to 20C0 m for temperature, water 
content, wind, and TBL thickness 
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For the days of the tests, 3-8 August 1987, profiles of the turbulent 
boundary layer thickness were obtained from the Blackadar model. One 
can see the thickness of the boundary layer generally grew during the 
day between 7:00 AM and 2:00 PM. Because of the meteorological 
conditions present, the boundary layer grew much more on 3 August 
and 4 August than it did on the other days of the tests. 


TURBULENT BOUNDARY LAYER THICKNESS 
ESTIMATION USING BLACKADAR MODEL 


Local time 

05:00 AM 07:00 AM 09:00 AM 11:00 AM 01:00 PM 03:00 PM 



300 400 500 600 700 800 900 

Time past midnight (minutes) 
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Before continuing, the analysis procedure will be summarized. On one 
track the BEAR sonic boom data was obtained and calibrated. From this 
data pressure versus time plots were obtained, which were subsequently 
sorted into waveshape categories. Given these plots, the sonic boom 
distortion quantification algorithm was run, and mean-squared 
deviations from the computer generated corresponding ideal N-waves 
were obtained for all the the waveforms. 

On the other track, the aircraft flight parameters and geometry were 
combined with the meteorological data and subsequent predictions of 
the boundary layer height from Blackadar's model. From this 
information the path length through the turbulence was found for 
each recorded sonic boom waveform. 

It was then possible to determine if a statistical correlation existed 
between the mean-squared deviation and the path length the sonic 
boom traveled through the turbulence. 


ANALYSIS PROCEDURE 
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During the analysis procedure it became immediately apparent that 
those waveforms which were either shaped as a Classic N-waves, or 
were nearly shaped as such, were primarily manifest only in the early 
morning hours. 

This result leads us to believe that as the turbulent boundary layer grew 
through the day, that the number of undistorted waveforms decreased. 
This result is averaging over all of the usable observation data. 


SIMPLE DEMONSTRATION: FRACTION OF RECORDED 
SONIC BOOM WAVEFORMS THAT ARE LEAST-DISTORTED 



445 recordings from 7/3 1 - 8/7 
overflights during 7 AM- 1PM 

Fraction is number of classic N 
and nearly classic N waveforms 
per number of waveforms 
recorded during the hour intervals 

classic N •+ least distorted 


number of 

recorded 

waveforms 


TBL thickness Increasing C> 
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This is a plot of the mean-squared deviation as a function of altitude (or 
roughly Mach number, since faster planes generally flew higher). It is 
apparent that increased altitude and speed imply decreased waveform 
distortion and a decreased spread of data points. Now higher and faster 
flying planes generally will have shorter propagation paths through 
the turbulent boundary layer, which can be shown from simple 
geometry. Thus, it appears as if longer propagation paths through the 
turbulence result in larger mean-squared deviations, i.e., more distorted 
waveforms. 


ALTITUDE AND MACH NUMBER INFLUENCE UPON 
WAVEFORM DISTORTION 


Increased altitude and speed ’*► decreased waveform distortion and 

decreased spread of data points 


F-18, F-15, SR-71, 
F-4, and F-16 data 
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Here are some results listed by plane type of the linear correlation 
coefficients between the mean-squared deviation and the path length 
through the turbulent boundary layer. It is seen that there is strong 
correlation in some cases ( F-18 and F-15) and fair correlation in the 
others. This also gives us some evidence that increasing distortion is 
correlated with increasing path length through the turbulence. 

The F-4 and F-16 data have the lowest correlations and the lowest 
altitudes. For these cases the sonic boom signatures probably did not 
have time to develop into an N-wave by the time it began to interact 
with the turbulence. The other aircraft types flown had too few data 
points to draw any statistical conclusions. 

We also are currently working on obtaining correlation coefficients 
grouped by altitude and mach number as well as by plane type, to 
determine how these factors interrelate. 


CORRELATION COEFFICIENTS FOR 
MEAN-SQUARED DEVIATION AND 
TBL PATH LENGTH 


Aircraft type 

Correlation coefficient: MSD vs 
TBL path length 

Number of samples 

F-18 

0.712 

62 

F-15 

0.591 

75 

SR- 71 

0.398 

48 

F-4 

0.324 

46 

F-16 

0.318 

65 


(Linear regression correlation coefficients) 


Strong correlation in some cases, fair correlation in others 


F-4 and F-16 have lowest correlations and the lowest altitudes 

- Sonic boom signature before TBL probably not N-wave 

- Less distance for nonlinear steepening to work before TBL 
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The conclusions of this study are the following: a strong linear 
correlation exists between the mean-squared deviation and path length 
through the turbulence for the F-18 and F-15 sonic boom data. Fair 
correlation exists for the SR-71, F-4, and F-16 data. An increase in 
altitude and speed results in decreased waveform distortion and a small 
deviation between distortion values. Looking at the waveform 
classification results, the large percentage of Classic N-wave data 
during each day’s early flights seems to correspond with the thin 
boundary layer at that time. 


Conclusions: 

• A strong linear correlation exists between mean-squared deviation 
and path length through the turbulence for F-18 and F-15 sonic 
boom data. 

• Fair correlation exists for SR-71, F-4, and F-16 data. 

• An increase in altitude and speed results in decreased waveform 
distortion and a smaller deviation between distortion values. 

• Looking at waveform classification results, the large percentage of 
classic N wave data during early flights seems to correspond with 

the thiii boundary layer at that time. 

pennSmte 
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The Implications of this study for high speed research are the 
following: Increased interaction between real atmospheric turbulence 
and actual sonic boom data does imply more distorted waveforms. The 
common assumption prevailing in the sonic boom propagation 
community for the last several years has been validated statistically. 
And most importantly, it is now clear that atmospheric turbulence will 
determine how well a shaped sonic boom will remain shaped as it 
propagates to the ground. We are now led to believe that higher and 
faster aircraft having shaped sonic booms will, on average, have more 
shaped boom preserved than will aircraft flying at lower altitudes and 
slower speeds, since flying higher and faster minimizes the path length 
through the turbulence. 


Implications for High Speed Research: 

• Increased interaction between real atmospheric turbulence and 
actual sonic booms does imply more distorted waveforms. 

• This common assumption has been validated statistically. 

• Most Importantly: Atmospheric turbulence primarily will 
determine how well a shaped sonic boom will remain shaped as 
it propagates to the ground. 


pennSkte 
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Abstract 


Theoretical research has been carried out to study the effect of free-stream 
turbulence on sonic boom pressure fields. A new transonic small-disturbance 
model to analyze the interactions of random disturbances with a weak shock 
has been developed. The model equation has an extended form of the clas- 
sic small-disturbance equation for unsteady transonic aerodynamics. An 
alternative approach shows that the pressure field may be described by an 
equation that has an extended form of the classic nonlinear acoustics equa- 
tion that describes the propagation of sound beams with narrow angular 
spectrum. The model shows that diffraction effects, nonlinear steepening 
effects, focusing and caustic effects and random induced vorticity fluctu- 
ations interact simultaneously to determine the development of the shock 
wave in space and time and the pressure field behind it. A finite-difference 
algorithm to solve the mixed-type elliptic-hyperbolic flows around the shock 
wave has also been developed. Numerical calculations of shock wave interac- 
tions with various deterministic and random fluctuations will be presented 
in a future report. 
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1 Introduction 


Experimental data exists showing the pressure profiles of sonic booms, cre- 
ated by the passage of a distant supersonic aircraft, can be drastically 
affected by free-stream atmospheric turbulence. 1-5 The usual N - wave or 
a shaped sonic boom profile can randomly exhibit either large pressure 
peaks with short rise times or rounded profiles with longer rise times or 
messy pressure signatures. Recent laboratory model experiments to study 
the effect of turbulence on the rise time and wave form of N - waves have 
shown similar results. 6 The interaction of the sonic boom with the atmo- 
spheric turbulence, specifically in the atmospheric boundary layer relatively 
near the ground, may result sometimes in higher, and may be unacceptable 
loudness levels.' Therefore, in order to get reasonable estimates of the sonic 
boom performance of various designs of a future supersonic transport air- 
plane it is essential to understand the basic interactions of the atmospheric 
turbulence with shock waves. 

The basic analysis of the distortion of sonic bangs by atmospheric tur- 
bulence was given by Crow. 8 Using a first-order acoustic scattering theory, 
Crow showed that the pressure perturbation behind the shock is related to 
the interaction of the shock with the disturbances it encounters while mov- 
ing in the atmosphere. The pressure profile can be calculated by a surface 
integral over a paraboloid of dependence, whose focus is the observation 
point and whose directix is the shock front. By describing the turbulent 
eddies in the Kolmogorov inertial subrange, it was found that the mean- 
square pressure perturbation behind the shock changes like (A p) 2 (t c /t) 7 ^ 6 
where (Ap) is the pressure jump across the shock, (t) is time after the shock 
passes an observation point and (< c ) is a critical time predicted in terms 
of meteorological conditions. Crow’s analysis predicts reasonable average 
values of the pressure fluctuations for times (f) comparable to (t c ). 

The singularity in the pressure perturbations near the shock front (when 
t — * 0) was analyzed by Plotkin and George 9 . A second-order acoustic 
scattering theory was used to describe shock rounded signatures. The av- 
erage of the diffraction effects was approximated a s a dissipation term. 
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The thickening of the shock is explained as a balance between nonlinear 
steepening effects and the dissipative effect of the turbulent scattering of 
acoustic energy out of the incident shock. Rise time predictions of this 
theory show some correlation with experimental data. On the other hand, 
Ffowcs Williams and Howe 10 examined the approaches that describe the 
possibility of a turbulent thickening of weak shock waves and reached a 
conclusion that atmospheric turbulence cannot be the cause of shock thick- 
ening. They suggested that weak shocks may attain a dispersed profile due 
to non-equilibrium gas effects. 

It should be emphasized here that the scattering analyses of References 
8-10 considered small turbulent perturbations against the shock strength, 
whereas in the case of the interaction of the sonic boom with atmospheric 
turbulence the flow random fluctuations may be of the same order of the 
shock weak strength and may strongly distort the shock front. The above 
analyses also did not account for shock jump conditions that must be sat- 
isfied in an inviscid analysis across any shock surface. The approximation 
made in Reference 9 of the average diffraction effects described as a dissi- 
pative term is also unclear. 

A different approach was taken by Pierce 11,12 . He interpreted the spikes 
observed on sonic boom pressure waveforms as being due to the simultane- 
ous focusing and diffraction of a nearly planar .V- wave by an inhomogeneous 
layer in the atmosphere. The shock front develops ripples that axe trans- 
formed into folds in the front when the shock passes vertices of caustics. 
This mechanism results in a fine structure of very small pressure jumps 
that correspond to the various segments of the folded wavefront. Pierce 12 
derived a stochastic model of a sharp shock propagating through a turbu- 
lent atmosphere to substantiate the very small discrete structure of sonic 
boom profiles. 

Pierce model 1112 , however, neglects nonlinear effects that become sig- 
nificant specifically near a caustic vertex as was shown by Cramer and 
Seebass 13 and Gill and Seebass 14 . Cramer and Seebass 13 described the fo- 
cusing of a very weak and slightly concave shock wave by the unsteady 
transonic small-disturbance flow equation. Gill and Seebass 14 derived an 
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approximate analytical solution of the steady transonic small-disturbance 
problem for the nonlinear behavior of a weak compression wave with a finite 
rise time near a caustic. They calculated the reflected shock wave from a 
caustic and provided an estimate of its strength. The experimental results 
of Sturtevant and Kulkarny 15 also show that focusing effects are specifically 
significant for weak shock waves as occurs in the case of the sonic boom 
signatures. 

Sparrow and Pierce 15 have recently presented a simple statistical predic- 
tion for how often sonic booms propagating in the earth’s boundary layer 
will encounter caustics. The theory is based on describing the variation 
of ray tube areas of a sound wave propagating in a turbulent medium by 
a generic harmonic oscillator equation. For realistic realizations of atmo- 
spheric turbulence the model predicts that sonic booms will exhibit spikes 
with the occurrence of caustics after a very short distance of travel in the 
random medium, thus agreeing with the predictions of Pierce 11 ' 12 . 

In a recent paper. Pierce 1, has derived a model equation to describe 
the development of sonic boom signatures in an atmospheric turbulence. 
The equation has been constructed by suing logical physical considerations 
only. It extends geometrical acoustic approximations to include convec- 
tion at the wave speed, diffraction effects, molecular relaxation, classical 
dissipation and nonlinear steepening effects. The atmospheric turbulence 
enters through an effective speed of sound which varies randomly in time 
and space. However, since this theory has not been developed consistently 
from the fluid dynamic equations, Pierce 17 raised questions whether all the 
effects axe necessary in his suggested model and how to accomplish a nu- 
merical or analytical solution to the problem. 

Related with the problem of the sonic boom interaction with turbulence 
is the basic question of the interaction of a shock wave with a vortical flow 
and specifically with a vortex or a train of eddies 18-21 . The research of 
the later problem was basically motivated by the interest to reduce the 
noise and vibrations produced by high-speed supersonic vehicles. For these 
problems, the interaction of relatively strong shock waves with turbulent 
jets or wakes is a significant source of noise. The shock-vortex system is a 
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basic element of these more complex interactions 21 . It can also shed light on 
the sonic boom interaction with atmospheric turbulence, specifically when 

the shocks are weak ana vortex strength is comparable with the shock 
strength. 


Experimental results of shock-vortex systems 18-21 revealed curved, dif- 
fracted shocks as well as complicate structures of reflected curved shock 
waves from the incident shock front due to the vortex induced flow field. 
The pressure field behind the shock contains regions of compressions and 
rarefactions that produce acoustic waves. Similar shock structures were 
also observed by Sturtevant and Kulkarny 15 who investigated the focusing 
of weak curved shock waves. Of specific interest are Dosanjh and Weeks 19 
measurements of the interaction of a shock wave with a vortex street. The 
shock front is distorted by the wake flow which probably results in a focusing 
process, while the vortex street is rapidly dissolved by the shock. 

The analyses of the shock-vortex interaction are limited to linear per- 
turbation theories only . 22-26 These analyses considered the jump conditions 
across a shock surface and predicted the development of vorticity waves, 
entropy waves and acoustic waves behind the shock front . 22-24 The acous- 
tic wave was approximated by a quadrupolar 26 or as a sum of monopole, 
dipole and quadrupole acoustic sources. However, since all of these theo- 
ries are linear, they cannot account for any nonlinear effects due to shock 
large distortions, focusing and caustic effects or nonlinear steepening ef- 
fects that are found in experiments 18-20 or in recent numerical simulations 
of shock- vortex interactions . 27-29 

The review of experimental and theoretical investigations of the interac- 
tion of shock waves with free stream vortical or turbulent flows shows that 
this complex nonlinear interaction is still an open problem . 27 Specifically, 
the improved simulation of sonic boom propagation through the real at- 
mosphere requires a better understanding of the interaction of sonic boom 
with atmospheric turbulence . 16 

The analysis of the experimental data and the theoretical approaches 
shows that in the case of the sonic boom, the shock waves near ground 
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are very weak, but still stronger than any acoustic wave. Also, flow fluc- 
tuations due the atmospheric turbulence can become comparable to the 
shock weak strength such that locally the shock strength can be strongly 
reduced or magnified and the shock wave front can be distorted significantly. 
Therefore, linearized acoustics and its second-order scattering problem, or 
first-order linear theories of shock- vorticity interaction do not represent cor- 
rectly the development of the weak shock and the pressure field behind it 
(see also Section 2). However, in a coordinate system moving with the basic 
weak shock, the problem may fit the transonic framework. 

This paper presents a new transonic small-disturbance model that has 
been developed to describe the interactions of random fluctuations with a 
weak shock wave. The model equation is found to have an extended form 
of the classic nonlinear acoustics equation that describes the propagation 
of sound beams with narrow angular spectrum (KKZ equation). 30,31 The 
model shows that diffraction effects, nonlinear steepening effects, focus- 
ing and caustic effects and random induced vorticity fluctuations interact 
simultaneously to determine the development of the shock wave in space 
and time and the pressure field behind it. A finite- difference algorithm to 
solve the mixed-type elliptic-hyperbolic flows around the shock wave is also 
presented. The results of the numerical calculations will be presented in a 
future report. It is expected to find solutions that will describe both peaked 
or rounded or messy pressure signatures as were recorded in experiments. 


2 The Breakdown of the Linearized Theory 

2.1 The Linearized Theory 


An inviscid and a non heat conducting flow is assumed. A normal shock 
with a uniform supersonic oncoming stream and a uniform subsonic outgo- 
ing flow is considered. The upstream flow ahead of he shock is characterized 
by a speed ( Uo a ) in the x- direction, pressure (poo) and density (p 0 a) and 
the downstream flow behind the shock by (t'oft). (pot) and (po*) respectively. 


71 


Assuming the shock front is given by the x = 0 plane, the jump conditions 
across a normal shock 32 show that: 


POaUoa = pObU( 


Ob 


POa b Oa + Poa — + Po6 


(1) 


,-yPOaU Q a + , POat- Oa — p06 f (}6 d* Pob^-'Ob 

- 7 - 1 1 7~1 

Small disturbances are considered in each of the uniform streams. The 
velocity vector (V’). pressure (P), density (p) and vorticity (w) are given 
ahead (j = a) and behind the shock (j =6) by: 

1 = £ o j( e + ev +•••) 

-j ~r ~lj 

P ; = Poj( 1 + ep\) + ■ • ■) (2) 

P) = Poj( 1 + £P\j + • ■ - ) 


<*-’ = eu; + • • ■ , w ,= \ 7 xu 

*" l; ~i^ 


Here ^^P\), P\j are functions of ( x,y,z,t ). An axial coordinate moving 

with the uniform speed is considered in each region, = x — U 0 j t,T = 

Uo } t. The substitution into the continuity, momentum and energy equations 

results to the leading order in ( V = 77“ e + ^-e + -f e ): 

'W ~x •■’y-y 


dr 


+ V 


v 


0 


<9u 


+ 


dr 7iVf 0 2 ; ~ 


VPu = 0 


( 3 ) 
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0 


dp\j 

dr 


+ 7(V-f ) = 

- ~i j 


Equations (3) result in: 


dpij dp\j 
dr 7 dr 


doj 



- V 2 Pi: = 0 


(4) 


where Mq } LqJc Qj and Cq ; — ipoj/poj- Equations (4) show that u> = 

~i/^’ 2 ) an< ^ ^at t ^ le pressure perturbation p Xj is described by the acous- 
tics equation. Therefore, the first order disturbance flow can be split into 
a linear combination of rotational and irrotational parts: v = v + v 

The rotational part can be described essentially by incompressible flow 
equations: 


Y'~V = 0 ' < 5a ) 

The irrotational (potential) part may be described by acoustics equations 
relative to the basic flow in each region: 


v 

~1;0 



Pi; = 


dr 


-V/, 


2 &<t>, 


Oj 


dr ^ 


- X7 2 <t>> = 0 


(5b) 


The first-order perturbation theory also considers the distortion of the 
shock front. Assuming that the perturbed shock front is given by x - 
£9i{y,z,t) = 0, the exact jump condition across the shock 32 result to the 
leading order in a set of conditions that must be satisfied along the x = 0 
plane for any (</, z, t): 


M0,y,z,O + ttu(0,y,j,O- “ = pu(0.y.:.t) + u lb {0,y,z,t) - 

U 0 b 


Uoa y, z,t) + 2u lo (0, y,s,t)) + y, z, t)j = 

= Uob ^Pu(0, y,z,t) + 2uii,(0. y, z,t) + — — p\t>{0, y, z, t)j 
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(p\b{0.y,z,t)+ 


+ u lt ( 0. y. --.<)) - (l + ~ ~T7TT2 
£ 06 \ ( 7 1 )-^*06 


^oa( ^ia(0,y,c,i) + yi y ) — 6 06(^16(0, j/, z, t ) 4- 9\ y ) 


Uo a {w Xa (0,y,:,t) + g u) = U 0 b{ u’u(0, y, 0 +yu) (6) 


Here y u = dg\/dt and (u tj , t’i ; , u,’i 7 ) axe the components of the velocity 
perturbation v . The linearized jump conditions in Eqs. (6) include the 


entropy increase produced by the shock. It can be shown from Eqs. (5) and 
(6) that, using the solution of the downstream equations for the disturbance 
flow, the shock conditions are adequate to describe the flow downstream and 
the disturbed motion of the shock wave for given upstream disturbances. 


2.2 One Dimensional Flow 


In the case of a one-dimensional flow, the rotational part vanishes identically 
and the solution of the acoustics equation ahead ( j = a) and behind ( j = b) 
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the shock is given by: 


= 4 -^) + 4 + 4 ) 

= Mo, (f -c(*, + ^) +,.,(*, )) 

where F and G are arbitrary functions that describe the upstream and 
downstream acoustic waves. Here p s is an arbitrary function that describes 
entropy waves that are convected with the flows (the entropy first-order 
disturbance is given by S ; - S 0j = es Xj = -ec p M 0jPsj (^), where 5, is 
the entropy in region j % 5 0<7 = c u ln(P 0; /po ; ) and c v ,c p are the specific 
heat constants). Assuming that no upstream acoustic waves can develop, 
specifically not in the flow behind the shock then G = 0. Then the shock 
jump conditions (6) provide a system of 3 linear equations for the solution 
of the downstream acoustic and entropy waves and and 

the shock position rate of change in time gu(t) in terms of the given acoustic 
and entropy perturbations F a (£ a - and p ja (£ a ) in the upstream flow. 
Let, 


FaO F a ( (Lqo “b Cg a )£), PsaQ — Paai 


( 8 ) 


FbO — Fb(~(Uob + Coi)#)- PsbO = Psb( “£^ 06 *) 


then from Eqs. (6) - (8) we get: 
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The determinant of the system (9) may be written in the form A = 
(Mq 0 — l)/n(Mo a ). Therefore, the solution of Eqs. (9) shows that the 
shock front motion and the perturbed flow behind it may be described by: 

* = vts — rgi(<;M 0a ) 

M) a - 1 


Uit = u 0i [ 1 + 


ML - 1 


u t (x, t: Moo) + 


Pn 


P\b 


Pot 1 + 


•v«. - 1 


p t (T,t;M 0a ) + 


) 


Pot 1 + 


ML - 1 


Pi(x, t: M 0a ) + 


(10) 


The functions §\,ui,p\,Pi can be expressed in terms of the given flow per- 
turbations ahead of the shock wave.. In principle, these expressions may 
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enable a spectral characterization of the pressure fluctuations and the tur- 
bulence downstream of the shock wave in terms of the spectral character- 
ization of the incoming turbulence. However, Eqs. (10) show that the 
linear approach is a non-uniform approach when the shock wave is very 
weak A/ Oo ► 1 + and the flow fluctuations are of the same order of the 
shock strength, e ^ — 1), as is the case of the interaction of the sonic 

boom with atmospheric turbulence. A similar nonuniformity is also ex- 
pected from the analysis of two- or three-dimensional flows. However, this 
uniformity problem leads to a different approach to study the interaction 
of a weak shock with comparable random fluctuations in the flow. 


3 A Transonic Small Disturbance Model 


The analysis of the linearized problem of the interaction of a weak shock 
with small disturbances shows that it is an invalid approach when the flow 
perturbations are of the order of the shock strength. Therefore, a different 
approach has been developed to study the interaction of weak shocks with 
comparable random fluctuations in the flow. In a coordinate system moving 
with a basic given weak shock, the problem may fit the transonic theory 
framework. A transonic small-disturbance model is developed to analyze 
the flow across a basic weak shock running in the ( — .r ) direction. A coordi- 
nate system attached to the basic shock is considered. The velocity vector 
(JO* pressure ( P ), density (p) and vorticitv (u. 1 ) are described every where 
in the flow by: 

V = Uoo{i(l + e 2/3 u + eui + e 4/3 u 2 + • ••) 

+>(«>i + e'* /3 L '2 H ) 4- k(ew { + e 4/3 u? 2 + •••)} 

P = Poo(l + e 2/3 p + epi + e' l/3 p 2 + • • •) (11) 

P = Poc( 1 + e 2/3 p + ep x + e 4/3 p 2 + • ■ • ) 
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^ j Wy + kiv z ) + c( 1 + J + k'y^zi ) 4* • * * 

** *■*’■''■' *v "w 

where U^, == a^,(l + f e 2 ' 3 ) is the speed of the basic shock (K > 0) and 
^OO 1 poo 1 Poo are the speed of sound, pressure and density of the unperturbed 
flow ahead of the shock. (e 2/3 ) represents the scale of strength of the basic 
weak shock where e 1. A rescaling of the x-coordinate and time (t) 
has also been considered: j* = and t’ = ta^e 1 ? 3 , such that each of 
the terms in (11) is a function of (x*,y, $,£*), The rescaling in x means 
a stretching of the picture of the flow around the basic shock in order to 
capture the basic nonlinear effects that occur in the flow across the shock. 
The rescaling in time accounts for low-frequency unsteady perturbations in 
the flow. The constant I\ reflects that the speed of the basic shock wave is a 
little higher than the speed of sound, ahead of the shock. The substitution 
of Eqs. (11) into the continuity, momentum and energy equations results 
to the leading orders in: 
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dx m dx" 

From the equation of state and the definition of entropy it can be shown 
that the temperature T m and entropy S are given by: 


= T^(l +t 2 ' 3 T + ■■■) 

S = S^( 1 + + ■ • • ) (17a) 

where 


7\o = P-JRpTC. 



c 

T = p - p, i = 7T L (p-7p) 
Equations (12) through (17) result in: 

u + P = f(y.:>t m ) 
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(13d) 


du dii'i 1 dg 

dz dx‘ 7 d: 

du dv\ _ 1 dg 

dy dx‘ 7 dy 


(18e) 


T = g-f-h- l)u, ( lSf ) 

where / and g are random induced fluctuations due to the free turbulence. 
The function g is related to the vorticity fluctuations in the flow. Equations 
(18) show that the axial perturbation (u), pressure perturbation (p) and 
density perturbation (p), that are of order of the shock strength (e 2 ^ 3 ), 
interact with the transverse velocity perturbations iq and uq. that are of a 
smaller scale (e). 


The substitution of u = g / 7 + u in (18c), (18d) and (18e) results in a 
problem for solving a velocity potential function o{x m ,y.:,t m ) where: 


u 



do 



do d<p 



(19a) 


2<£ r *C -+* A + b / + ( 7 "t" 1 )©r* I <5 r *x‘ — ( Oyy + <t) ZZ ) — 
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In a conservative form Eq. ( 19b) is given by: 


^2<^ r . + — + (( A + g/“t + f)O x • + ( 7 + 1 )<2>j«/2))x* — (^yly — (<^r)x — 0 

(19c) 

The exact shock jump conditions (Ref. 32) must be satisfied along any 
shock surface x’ — h{y,:,t') = 0 that may appear in the solution. To the 
leading orders they result in: 

[/] = 0. [g] = Q (20a) 
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,dh 


r i On 

_ “[^ x ‘] dt* + (^" + 5 F + /H <i, r‘] + (7 + l) 


. .dh , .dh 
+ [M^+(^]^r = 0 (20b) 


[0y] "b (0r *] Qy ~~ ^ [0*] "b [0x*]^j = 0. (20c) 

where [a] represents the jump across the shock property a, [a] = a g — a A . 
Equations (20a) show that to the leading order there is no jump in entropy 
across the shock, [5] = 0. Equations (11) and (IS) also show that the local 
Mach number Mr at any point in the flow is given by: 


M} - 1 = e 


2/3, 


11 — { (7 + 1 )Ox* # + A + / + 


7} 


( 21 ) 


The flow is locally supersonic when (7 + l)0 r . -f A' -f / 4- * >0, sonic when 
(7 + 1)0*. + A' + / + * = 0. and subsonic when (7 4 - 1)0*. + A' 4- / + * < 0. 
Equations (19) and (20) are an extended version of the classic unsteady 
small-disturbance equation for transonic aerodynamics (Cole and Cook 33 ). 
The only changes are due to the random terms g and /. Starting from 
given functions for / and g and initial conditions that describe a given 
basic shock, Eqs. (19) and (20) can be integrated in space and time to 
describe the development of the shock wave and pressure field behind it. A 
numerical algorithm to solve these equations is described in Section 4. 


An alternative approach may be found by taking a x'- derivative of ( 18c) 
and using Eqs. (18a) and (14). The pressure perturbation (p) satisfies the 
equation: 


d ( dp K + / -t- gj_ 7 _ 7 + 1 dp \ _ 1 / cPp d 2 p \ 
dx‘ \df 2 27 P dx’ ) ~ 2 \dy 2 + dz 2 ) 


Equation (22) is an extended version of the classic KKZ equation that 
describes the propagation of nonlinear sound beams with narrow angular 
spectrum in an inviscid fluid (Zabolotskaya et al. 30 . Kuznetsov 31 ). Eq. 
(22) also has a similar form to the model equation that has been recently 
developed by Pierce 1 ' using logical considerations only. 
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Equations (19) and (22) show that diffraction effects, nonlinear steep- 
ening, focusing and caustic effects and random induced fluctuations due 
to turbulence interact simultaneously to determine the development of the 
shock wave in space and time and the pressure field behind it. Turbulence 
tends to change the local speed of sound in the flow across the shock and 
through this effect to reduce or to magnify the strength of the jump along 
the basic shock (see Eq. (21)) or to distort the shock front. These changes 
may result in unsteady motion of the shock front or in caustic vertices or 
in reflected shocks behind the incident wave that can produce the variety 
of pressure signatures of sonic booms that are measured in experiments. 


4 Finite Difference Scheme 


A finite difference algorithm to solve the unsteady mixed-type elliptic- 
hyperbolic flow around the shock wave has been developed. Murman and 
Cole 34 and Cole and Cook 33 techniques are used. A fully conservative 
scheme that is based on the conservative form (Eq. 19c) is used. In this 
way the difference equations also contain the shock relations (Eqs. (20)). 


Consider a uniform finite difference mesh ( Ax'. Ay, Ax, At’) in space 
and time, with points (x m ,y,;,t m ) labeled by ( i.j.k,n ). The results can be 
easily generalized to a variable mesh. Eq. (19c) can be expressed in a con- 
servative flux form for a box centered on a mesh point (i, j, k). Therefore, 



+ 


Ax* 


K +-g + f ) 0 X . + (t- + 1)<P 
t 7 / 


1 


^.-j J 


— ^A + —g + J J + (7 + 1)0 

~ )(i.j - i.fc.n) | — ^t|(^- ^.n) ~( < Pz )( t.j.fc- 1 ,n) j ~ 0 

(23) 
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{4>y) a-nd (<£,) are always calculated from a centered expression. However, 
the approximation of [<p z ) strongly depends on whether locally, at a point, 
the flow is subsonic, supersonic, sonic or it is a shock point. Extending 
References 33, 34 methodologies to our case and using Eq. (21), a centered 
approximation and a backward expression are given for u': 

u U.j,k,n) = & + f(j,k,n) + -g(j, k.n) 

7 

7 + 1 

+ o A r - + n ) ~ o(i - l.j. k,n)) 


u U.j,k,n) - h + f(j.k.n) + -g(j,k,n) 

7 


7 + 1 


+ - <p(i - 2 ,j,k, n)) 


(24) 


The local type of the flow is determined by the following table: 33,3,4 


condition 

u mc 

u mb 

local flow is 

1 

< 0 

< 0 

subsonic 

2 

> 0 

> 0 

supersonic 

3 

> 0 

< 0 

a sonic point 

4 

< 0 

>0 

a shock point 


Table 1 


Eq. (23) is developed in a specific form according to the local type of the 
flow. When the flow is locally subsonic, an elliptic difference form is used: 


83 



G(i,j, k. n) 


- — ~-(0(t + 1 ,j,k.n) ~ 0(1 - 1 .j,k,n - 1)) 
Ax* A t m 


2 1 

— n - 1) + -£j;( 9 (j,k,n) - g(j,k,n - 1)) 


+ (/\ + f{j. k. n) 4- ~g(j. k. n ) 


+ 


T~r— {0(i + l,j, k. n) - <?{i - l,j, t, n)) > 
2 A x' J 


o(i 4- 1 ,j,k,n) — 2o( i,j,k.n) + d>( i — 1 ,j,k,n) 
(Ax-) 2 


(Ay) 


~(o(i,j + l.k\n)- 2 k.n) + o(i.j - l,k,n)) 


(Ax) 2 


{<p(i.j,k 4- l,n) — 2 ®{i.j<k<n) 4- <p(i,j,k - l,n)) = 0. 


(25) 


When the flow is locally supersonic, a hyperbolic difference form is used: 

2 

G(i,j,k,n) = " ( <?(i.j . k, n ) - 6( i - 1 . j, k. n ) ) - 


o 1 

~0 z .(i.j,k, n - 1) + ^ - {g(i>J,k,n) - g{t,j,k,n - 1)) 
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+ { K + f(j , k.n)+ - g(j , fc, a) 


7 + 1 

n) - o(i - 2 J,k.n)) 


Ar, n) - 2d>(z - 1 ,j, k, n) + 4 »(t - 2, y, fc, n) 

(Ax’) 2 


(Ay) 


■j( <z>( i. y + 1 . k. n) — 2o (i.j. k, n) + <p(i. j — 1 .k,n)) 


(Ar) 


~(q{i. j. k + 1 . n) - 2 6(i,j, k, n) + <p{i,j, k - 1 , n)) = 0 . 


When the flow is locally sonic. (7 + 1 )<z> r . + I\ + / + = 0. Then the sonic 

point difference form is: 

G(i,j, k,n) = * — (0(ij + 1 .k. n ) - 2 o(f.y, k. n) + o(i,j - 1 . k, n)) 

(Ay) 2 


+ (A^ i0{l ' J ' k + ~ + <P{i>j,k - 1 , n)) 


(7 + 1 ) A t‘ 


7 - 1 


(gij.k.n) - g(j.k\n - 1 )) 


-2(/(y\ k.n) A\ n - 1))J =0. (27) 

When locally there is a shock point, a shock point difference operator is 
used where the flux ahead of the shock may be approximated by a backward 
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formula and the flux behind the shock by a centered formula. 


G(i, j. k, n) 


Ax-A r (o(/ + lj - k '- n) ~ ~ W- k - n )) 

2 1 

~AF 0l ' (i ' ;,t ' n " U + 7^7 z(9U'k<n)-g(j,k,n- 1)) 


+ | A' + f(j, k\n) 4- A\n) 


7 + 1 

+ ^ — ( °i i + 1. J* k. n ) - o(i.j.k.n) + <p(i — 1 ,j,k.n) 


+<->( i - 2.j, k. n ))| • 

<p(; + 1 J,k.n) - o(i,j.k\n) - <j>{i - l.j, fc. n) + -2 J,k,n) 

(Ax *) 2 


1 

~ ( Ay)^ 0 ^' 7 + 1 '^' n ^ - n ) + $>(*J ~ 1, fc, n)) - 

1 

-^(CK,;. fc + l.n) - 2o(t.j. k. n) + A: - 1, n)) = 0. 

(28) 


Starting from initial conditions that describe a given shock wave in the 
space for t = 0 (or n = 0), and given the functions f{y,zA) and g(y,z,t), 
equations (25) through (2S) can be applied for n = 1 at any mesh point 
according to Table 1. They can be solved by an iterative point or line - 
or plane - relaxation algorithm until at any point max |G( t, j, k, 1 )| < 6 
where 6 is a given small tolerance of convergence. Then o r .(i, j, k, 1) can 
be calculated at any mesh point and the process is restarted for the next 
time step. In this way the shock motion and pressure field behind it can 
be integrated in space and time and the effect of various deterministic and 
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random fluctuations / and g can be studied. Numerical calculations of 
various examples are underway and will be presented in a future report. 


5 Conclusions 


The review of the theoretical studies of the interaction of shock waves with 
free stream vortical flows or turbulence shows that this complex nonlinear 
interaction is still an open problem to analyze. The analysis of the lin- 
earized problem of the interaction of a weak shock with relatively small 
disturbances shows that it is an invalid approach when the perturbations 
are of the order of shock strength. However, in a coordinate system mov- 
ing with the basic weak shock, the problem may fit the transonic theory 
framework. 

A new transonic small-disturbance model has been developed where a 
rescaling of the axial coordinate and time has been considered to capture the 
basic nonlinear effects that occur in the flow across the shock. This model 
results in two alternative approaches: (1) an equation for solving a velocity 
potential function that is described by an extended version of the classic 
small-disturbance equation for unsteady transonic aerodynamics; 33 and (2) 
a nonlinear equation to describe the pressure field that is similar to the 
model equation recently presented by Pierce 1, using logical considerations 
only. This equation also has extended form of the classic equation that 
describes the propagation of nonlinear sound beams with narrow angular 
spectrum. 30 31 

Both approaches show that diffraction effects, nonlinear steepening, fo- 
cusing and caustic effects and random induced turbulence fluctuations in- 
teract simultaneously to determine the development of the shock wave in 
space and time and the pressure field behind it. Turbulence fluctuations 
tend to change the local speed of sound in the flow across the shock and 
through this effect to reduce or magnify the strength of the basic shock. 

A finite difference scheme that uses Murman and Cole 34 finite-difference 
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techniques for solving mixed- type elliptic-hyperbolic flows with shock waves 
has also been presented. Numerical calculations of the interaction of shock 
waves with various deterministic and random fluctuations will be presented 
in a future report. We will also look for analytical methods to identify the 
basic relations between given turbulence properties and the development 
of the shock waves and the pressure field behind it. It is expected to 
find solutions that will describe both peaked or rounded or messy pressure 
signatures as were recorded in experiments. We also intend to extend the 
model to include humidity and winds effects that axe also known to have a 
significant effect on sonic boom pressure peaks and rise times. 
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In this work, a ray theory approach is used to examine the propagation of sonic 
booms through a turbulent ground layer, and to make predictions about the 
received waveform. The rays are not propagated one at a time, as is typical in ray 
theory; instead, sufficient rays to represent a continuous wave front are 
propagated together. New rays are interpolated as needed to maintain the 
continuity of the wave front. In order to predict the received boom signature, the 
wave front is searched for eigenrays after it has propagated to the receiver. 


OVERVIEW 

• Rays describing a wave front propagate 
through an Instantaneous "snapshot" of the 
turbulence. 

• Turbulence produces focusing and defocusing 
of portions of the wave front, which results in 
caustic formation, wave front "folding", and 
multiple eigenray paths to the receiver. 

• The eigenrays to the receiver are identified. 

• The respective arrival times and ray tube 
areas of the eigenrays, along with the 
identification of caustics, generate 

the predicted waveform at the receiver. 

• If repeated many times, this generates a 
statistical description of the predicted wave 
form characteristics. 
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The Comte-Bellot turbulence model (Ref. 1) is used to generate an instantaneous 
"snapshot" of the turbulent field. The transient acoustic wave is assumed to be 
sufficiently short in duration such that the time-dependance of the turbulent field 
may be neglected. 


Turbulence Model 
(Comte-Bellot ’91) 

Instantaneous realization of incompressible, 
isotropic turbulence is represented by a sum of 
Fourier modes: 

N 

w(x) = X aj cos(kj-x + <|>j) 

j = i 

where directions of aj and kj are random with 
the provision that 

(aj-kj) = 0 for each j 

For a given mode kj, the magnitude laji is given 
by 

lajl ~ >/E(k) 6k 

5k is the separation between modes. 

The spectral energy density E(k) is given by the 
Von K£rm£n model: 

E(k) - ^exp(-2.25(nk)4/3) 

( k2 + i ) 17 ' 6 

Lo = integral length scale 
rj = Kolmogorov scale 

The magnitude of the rms velocity is related by 

Ivrmsl = \ Zlal 2 

j = i 
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In this model, we use60 Fourier modes logarithmically distributed between 
wavenumbers 10 2 and 10 (nr 1 )- The integral length scale was chosen to be 
100m, as was the thickness of the turbulent layer. The rms wind velocity was 
chosen to be 1 m/s. This corresponds to a mild turbulent layer, such as might be 
found in the morning on a clear day. These values are used for all of the 
remaining figures and discussions. 


Turbulent layer Energy Spectrum 



Lq = 100 m r] = 0.01 m 

v™ = 1 m/s 

Used 60 modes, from k = 0.01 to 10 m 1 
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A number of rays, with the same starting conditions, are propagated through 
different realizations or "snapshots" of the turbulent field. Each ray will be 
displaced by the turbulence away from the undistorted ray path, which in this case 
would be represented by a horizontal line. 


Ray Paths 



• Each ray represents a different realization of the 
turbulent layer. 

• The turbulence displaces each ray from the 
horizontal (undistorted) ray path. 


97 



If a large number of rays are propagated through different turbulence realizations, 
a pattern or distribution of ray displacement may be developed. In our case, 95% 
of the rays fall within a circle of radius 2 meters around the zero-turbulence ray 
path. By symmetry arguments, this means that any eigenrays have a 95% 
probability of starting within 2 meters of the zero-turbulence eigenray. This 
statistical approach allows us to drastically reduce our eigenray search area to a 
feasible quantity. 


Ray displacement distribution after 100 m propagation 



X 


• 95% of the rays fall within a circle of radius 2 m 
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The development of the ray tube area along the ray paths gives an indication of 
how much the turbulent field is distorting the wave fronts. A ray crossing the 
horizontal as indicates that the ray has passed through a caustic at that point. 



Z 


• Each ray represents a different realization of the 
turbulent layer. 

• Crossing the horizontal axis indicates that the ray 
has encountered a caustic. 
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In this figure, a linear "slice" of a wavefront is propagated through 100 m of 
turbulence. In this case, the distortion is slight and no caustics are observed. 
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Again, a linear "slice" of the wavefront is propagated through a realization of 
turbulence. Although the statistical parameters are unchanged, in this case, the 
turbulence has a marked effect on the wavefronts. After 20 m, caustics begin to 
form which eventually overlap, producing, in the end, a highly folded wave front, 
with multiple eigenrays to the receiver. 
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It must b e re membered that the wavefront distortion is three-dimensional. The 
two plots below show the distortion of the original wavefront "slice" in the 
previous two figures, show in the plane normal to the direction of propagation. 
Note that the second figure shows considerable distortion due to the presence of 
numerous caustics. 
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For four realizations of the turbulent field, the wave fronts were propagated, 
eigenrays were found, and the resultant waveforms were calculated. The initial 
waveform was generated by the ZEPHYRUS model (Ref. 2); it represents a 
typical sonic boom waveform in the absence of turbulence. The next two plots 
demonstrate theresulting waveforms when the wave front is spreading, or 
defocusing, and when the wave front is focusing, but has not formed a caustic. 
The last plot displays the U wave resulting from multiple eigenrays, some of 
which have passed through one or more caustics. 



In nonlinear geometric acoustics, the effects of self-refraction may usually be 
ignored. Although nonlinear effects may displace a ray from the small-signal ray 
path, the properties of the wave front is usually slowly varying in the plane of the 
wave front and so, the equivalent nonlinear wave front is virtually identical to the 
original. As we've seen, however, in the case of propagation through turbulence, 
we've seen that the wave fronts may become very distorted and so the assumptions 
that lead to neglecting self-refraction must be examined more closely. This is 
most easily tested by comparing the same ray with and without the nonlinear 
correction. We first start with the nonlinear ray path equations given below. 


Nonlinear Ray Equations 

The ray path equations may be modified to include self- 
refraction as follows: 


dx 

dt 


( P'a\z 

[ Co+ Co*P 


(W + u’) + 


dp _ 

dt " 


f 

I- P,V'(W, + u'i) 

LV i 


(1 -W-p-u'-p)J 
P’a 


( Paw Pa\ 

l^vrhv)i" ; 

(1 - W p - u p) 


where 
a = 


(P-l)cp 
Po 

p is the slowness vector 

P' and u’ are the acoustical overpressure and 


and 


particle velocity 
V' is 


is the spatial operator in retarded time 
coordinates: 

V' = V + p| 
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The simplest wave front property to calculate is the ray tube area. When we 
compare the results for a number of different turbulence realizations, we see that 
the ray tubes with and without the nonlinear correction give almost the same 
result. This indicates that, for these environments, the wave fronts remain 
sufficiently smooth that we may continue to ignore the effects of self-refraction. 



• Initial 200 Pa acoustic overpressure 

• The nonlinear correction to the ray paths 
makes little difference. 
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The ray theory approach has been demonstrated to be a useful tool for the 
investigation of propagation through turbulence. The next step will likely be to 
attempt prediction for more severe turbulence, to see if waveforms of more 
complex structure that have been observed, such as multiply peaked or rounded, 
can be simulated by this method. 

It is fortunate that the nonlinear distortion of the ray paths may be neglected, as 
this simplifies the goal of sonic boom prediction. 


Conclusions 


• A ray theory approach provides a useful 
tool for investigating the properties of 
propagation through turbulence. 

• Wavefront folding and multiple eigenrays 
are good candidates for explaining some 
of the structure commonly observed 

in sonic booms. 

• Nonlinear distortion of the ray paths may 
be safely ignored. 
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A model experiment was designed and built to simulate the propagation 
of sonic booms through atmospheric turbulence. The setup of the model 
experiment is described briefly. Measurements of the N waves after they 
propagated across the turbulent velocity field reveal the same waveform 
distortion and change in rise time as for sonic booms. 

The data from the model experiment is used to test sonic boom models. 
Some models yield predictions for the waveform distortion, while others 
give estimates of the rise time of the sonic booms. 

A new theoretical model for the propagation of plane N waves through 
a turbulent medium is described. 


Introduction. 


- model experiment: - successful in simulating the propagation of sonic 

booms through atmospheric turbulence 
- setup and results 


- model experiment data is used to test sonic boom models 

- waveform distortion models 

- rise time prediction models 

- new theoretical model for the propagation of plane N waves through 
a turbulent medium 
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In the model experiment the N waves are generated by a spark source. 
The spark-produced N wave is a spherically spreading wave, but it is also 
possible to create a locally plane N wave by inserting a paraboloidal reflec- 
tor. The mirror is positioned so that the spark gap is at the focus of the 
paraboloidal reflector. 

A plane jet generates the turbulent velocity field. A centrifugal fan 
blows air into a plenum chamber. The jet is formed when the air exits 
the chamber through the nozzle. The jet nozzle velocity is controlled by a 
Variac variable voltage controller and by adjusting the width of the nozzle. 
The plane jet characteristics are measured by hot-wire anemometry. 

The N waves are measured by a wide band condensor microphone. Rise 
times as small as 0.45 //s can be measured. 
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Examples of waveform distortion are shown. The upper left signature 
is that of a reference plane N wave recorded in the absence of turbulence. 
All other signatures represent waveforms measured after the plane N waves 
propagated through the turbulent velocity field. The distortion of the wave- 
forms is similar as observed for sonic boom signatures. The distortion of 
the wavefront is most pronounced near the front and tail shocks. The fact 
that the distortion of the tail shock has the same pattern as that of the 
front shock is an indication that the turbulence is frozen during passage of 
the N wave. Variations in waveform from peaked to rounded and U-shaped 
are apparent. Double-peaked and multiple-peaked waveforms and messy 
wave shapes are also represented. 


Waveform examples 

Reference N wave 






The first sonic boom model we review is Crow’s waveform distortion 
model (Ref. 1). Crow’s model is based on first order scattering theory. He 
modeled the sonic boom as a step shock of strength A p. The mean-squared 
pressure perturbation equals (p|) = A p 2 (t c /t) 7 ^ 6 , where t c is a critical time 
that is a function of the turbulence characteristics. The graph presents an 
example of the variance of an N wave for a value of t c = 2 ms. A finite, 
very large value is obtained for the mean-squared pressure perturbation 
near the shock. Since the theory is a first-order scattering theory, both the 
incident and the scattered wave propagate at the ambient speed of sound 
Co- However, from geometric acoustics we know that some ray paths might 
exist along which the actual propagation speed is faster than the ambient 
speed of sound. If we want to compare Crow’s prediction with experimental 
data, then we have to shift the time origin of each sample waveform so that 
it begins at the time of shock arrival. 


Testing of previous models 
1 . Crow's distortion model 


- step shock :(^ 2 (0) = (tc/t) 7/6 t ) = p s /Ap, 

t c = - (sec 0) 1 1/7 [ I*' 1.33 Co~ 2 h 5 / 6 e^/ 3 dh 
Cq [Jo 

- Result of rms scattered pressure (tc = 2 ms) 



1 6/7 


- Problems: < V> 2 (0) > 10 8 

- incident wave of form f(x -f co t) 

- shock arrival time? 

- comparison with experiment? 
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Plotkin (Ref. 2) showed that an incident wave of arbitrary structure can 
be represented as a sum of infinitesimal step shocks. When the incident 
wave is modeled as a ramp shock instead of a step shock, an upper bound 
for the maximum pressure perturbation can be found, which is given by 

(p|) ,/2 = yAp , where to is the rise time of the ramp shock. 

Lipkens (Ref. 3) extended this model for an N wave. The rms pressure 
perturbation is presented in the first graph for an N wave with a rise time 
of 1 /zs and a critical time of 0.33 /is. In order to calculate the rms pres- 
sure perturbation for the measured waveforms of the model experiment, we 
shifted the time origin of each waveform so that the times corresponding 
to 50 % of peak pressure all coincide. 

A comparison between the measured distortion and Crow’s prediction is 
presented in the lower graph. The measured distortion has the same general 
behavior as Crow’s prediction, but the maximum pressure perturbation 
according to Crow’s prediction is larger than the measured one by a factor 
of more than ten. 


- Plotkin (1971): extension of Crow’s model for a ramp shock 

- incident wave of arbitrary structure: sum of infinitesimal steps 

- upper bound for rms pressure perturbation 




1/2 

max 



7/12 


- extension to N wave 



Time (ps) 

- measurements: shift each waveform — > 50 % peak pressure times coincide 



Time (ps) 


114 





It is now widely accepted that molecular relaxation, especially that of 
nitrogen, is responsible for the large rise time of sonic booms. Controversy 
still exists whether turbulence has a pronounced effect on rise time. When 
a theory for the propagation of shock waves and transients through a tur- 
bulent medium is developed, it is important to incorporate the effect of 
shock arrival time correctly. The rise time of the stochastic mean of a set 
of waveforms represents an insignificant upper bound to the average of the 
rise times of each individual waveform. The graph shows a simple example 
that demonstrates that the rise time of the stochastic mean waveform of 
five step shocks, each having a rise time of 1 /is, is more than tenfold the 
average of the rise time of the individual realizations. 


Rise time prediction models 

- Turbulence cause of large rise time? What mechanism? 

- Question of shock arrival time 
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Pierce’s model (Ref. 4) is based on the mechanism of wavefront folding at 
a caustic. If at a certain instant turbulence causes a ripple to develop on the 
shock front, a caustic is formed when the wavefront is propagated according 
to geometric acoustics. Inside the caustic three segments of the shock front 
arrive instead of one. Pierce argued that this process could occur many 
times if the turbulence intensity is large. A receiver then “sees” many 
segments of a multifolded wavefront at slighly different arrival times. The 
result is that instead of a sharp shock front a shock is received that consists 
of many smaller shocks at different arrival times. The overall result is a 
rounded shock front. The lower graph shows the mean waveform calculated 
according to Pierce’s theory. Again, t c is Crow’s critical time. 
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The mean waveform is expressed as a function of Crow’s critical time t c . 
The parameters E and bi are dependent on the structure of the atmospheric 
turbulence. 

Again, the individual waveforms of the model experiment are shifted 
in order that times corresponding to 50 % of peak pressure coincide. The 
plot compares Pierce’s prediction for the mean waveform and the computed 
mean of the shifted individual waveforms. A good agreement is reached. 

In order to confirm this correlation, we performed an experiment at five 
different jet nozzle velocities. The comparison between Pierce’s prediction 
and the measurements is fairly accurate. A maximum discrepancy of about 
30 % is observed for a nozzle velocity of 31.3 m/s. 


- mean waveform 

= exp {-Et 1/6 e- bli7/e } t = t/tc 
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jet nozzle 
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TPIerce 

fas) 

Tmeas. 
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12.4 

0.23 
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0.685 
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0.32 
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Kulkarny and White (Ref. 5) and White (Ref. 6) developed a model 
for the plane wave propagation through a 2-D and 3-D, random, isotropic 
medium. The model is based on geometric acoustics. The results shown 
here are for a 3-D, random medium. White derived a uniform probability 
density function p(t) for the occurence of a first caustic. The parameter t is 
a nondimensional variable. The only information needed about the random 
medium is its correlation function. Once this information is known, the 
scaling variable 7 is calculated. The graph shows the probability density 
curve. It is observed that the most likely position for the occurence of a 
first caustic is at t = 1.3. In the table, values for the most likely position 
of a first caustic and the mean distance to a first caustic are shown. As is 
noticed, it is possible that an N wave will pass through a caustic. However, 
it is unlikely that the wave will pass through more than one caustic. 


White and Kulkarny: plane wave propagation through a 3-D, 
isotropic, random medium (geometric acoustics) 

- probability density function p(t) for the occurence of a first caustic 

m = i = ■?!>., 



jet nozzle 
velocity 
(m/s) 

most likely pos. 
caustic 
m 

mean distance 
first caustic 
m 

12.4 

1.038 


18.3 

0.798 


22.7 

0.694 

0.918 

26.6 

0.619 

0.818 f 

31.3 

0.556 

0.735 | 
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Plotkin and George (Ref. 7) developed a model based on second order 
scattering theory. They derived a Burgers equation in which the absorption 
term is a function of the turbulence characteristics. L 0 is the integral 
length scale, a measure of the eddies of permanent character, and e T = 

^ (Ac + Au||) ^ /c^ is an effective turbulence Mach number. The rise time 

is determined by the balance between nonlinear steepening and scattering 
by turbulence. An expression for the rise time is obtained. Plotkin (Ref. 2) 
compared results from their model with measurements and obtained a good 
correlation. It is, however, difficult to obtain accurate estimates of the 
integral length scale and the turbulence Mach number of the atmosphere. 
A controversy still exists as to whether travel time variations are accounted 
for correctly or not. 
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In order to compare Plotkin’s prediction with measurements, we com- 
puted the mean waveform in two different ways. First, we computed the 
mean waveform by shifting the time origins in a similar way as described 
before and we call this the time shifted mean. Second, we calculated the 
stochastic mean waveform without any correction for the arrival time of 
each individual waveform and we call this the stochastic mean. The graph 
shows an example of the difference for the two computed mean waveforms. 
As is seen, the stochastic mean presents an insignificant upper bound to 
the rise time of the time shifted mean. 

Again we compare results for five different nozzle velocities. It is no- 
ticed that Plotkin and George’s prediction has a better correlation with the 
stochastic mean than with the time shifted mean. A conclusion that seems 
apparent from the results is that travel time variations of individual waves 
are not accounted for correctly in Plotkin and George’s model. 
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The last model is that of Ffowcs Williams and Howe (Ref. 8) (FfW & 
H). In their paper FfW & H mention that the Burgers equation derived 
by Plotkin and George represents the stochastic mean. FfW & H warn 
about possible misinterpretation of the results of the Burgers equation as 
an energy loss, while in reality it describes the loss of coherence of the 
mean wave because of the random convection of the shock fronts. The 
model by FfW & H is based on a multiple scattering theory. A diffusion 
equation is obtained that describes the acoustic energy £ in wavenumber 
space as a function of the turbulence Mach number m and a length scale A 
related to the Taylor microscale. An expression for the shock thickness 6 is 
derived as a function of the incident shock thickness <5 0 and the integrated 
scattering diffusivity fi. FfW & H found at most an increase of 30 % in 
the rise time and concluded in their paper that molecular relaxation must 
be the cause of the large rise times of booms. Plotkin (Ref. 9) argued that 
since his model does not yield an acoustic energy loss but just a spatial 
relocation, one would not expect a change in rise time according to the 


definition employed by FfW & H. 


3. Ffowcs Williams and Howe (1973): Multiple scattering 

- Plotkin’s approach describes stochastic mean properties of boom 

- Multiple scattering theory: diffusion equation for distribution of 
acoustic energy in wavenumber space 


d£{V) d£{V) Corn 2 k 2 2 C n ^ 

“9T + <: "-air = _ 2A- v ^ (k) 




<5 = 6 0 [1 + /i(x)] 

- At most a 30 % increase in rise time 

- Molecular relaxation is responsible for large rise time 

- Problem: phase scrambling is not accounted for 
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The upper graph presents a comparison of the energy spectrum of an 
ideal N wave, its rise time and duration equal that of the spark-produced N 
wave, with the averaged spectrum of 200 N waves recorded in the absence 
of turbulence. The energy spectrum of the spark-produced N waves closely 
resembles that of an ideal N wave. The troughs and peaks of the measured 
spectrum are more rounded than that of the ideal N wave. 

The middle graph shows a comparison of the averaged energy spectrum 
of 200 N waves recorded in the absence of turbulence with that of 200 N 
waves measured after propagation through the plane jet turbulence. Again, 
both spectra are very similar, and troughs and peaks are more rounded for 
the N waves that propagated through the turbulent medium. However, no 
significant redistribution of acoustic energy is observed in the spectrum, as 
was predicted by Ffowcs Williams and Howe. 

The table presents a comparison between the prediction of Ffowcs Williams 
and Howe and the measured values for the rise time. Ffowcs Williams and 
Howe’s model clearly yields values for the rise time that are much smaller 
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A new model has been developed for plane wave propagation through a 
statistically random, isotropic medium (Ref. 3). The random medium con- 
sists of a turbulent velocity field. A linear acoustic wave equation (Ref. 2) 
is derived in which first and second order turbulence effects are included. A 
perturbation scheme is used to solve the wave equation up to second order. 

The turbulent velocity field model was developed by Karweit et a 1. 
(Ref. 10). Von Karman’s model for incompressible, isotropic turbulence 
is used to obtain an expression for the 3-D turbulence energy density spec- 
trum. The spectrum is characterized by two length scales. L 0 is an outer 
length scale, and t? is the Kolmogorov microscale. 

The 1-D energy spectrum of the plane jet was measured by hot-wire 
anemometry. If we assume the turbulence is isotropic, the 3-D energy 
spectrum can be derived. A good agreement is reached between the model 
and the measurement. 


New numerical model for plane wave propagation through 
a statistically random, Isotropic medium 

- Turbulent velocity field model 

- Linear acoustic wave equation, second order turb. effects 

- Perturbation solution 

- Results 


1. Turbulent velocity field model (Blanc-Benon, Comte-Bellot, 1991) 


- Von Kdrmdn’s model for incompressible, isotropic turbulence 


E(*) = 


55 T(5/6) u ' 2 k 4 

9 VW (1/3) i^ + V) " /6 eXP 


(-2.25(7?*:) 4/3 ) 
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The turbulent velocity field consists of a sum of discrete Fourier velocity 
modes that arc randomly oriented in space. The wave vector geometry of a 
single Fourier velocity mode is shown in the first graph. The angles 6 and 
<j> determine the orientation of the wavevector k. The probability density 
function of both angles axe chosen in order to ensure statistical entropy 
with respect to k. With each turbulence wave vector, a velocity vector 
a ( k ) is associated. Because the turbulent velocity field is incompressible, 
the velocity vector lies in a plane perpendicular to k. The random angle tp 
determines the direction of a (k). The amplitude of a (k) is defined by von 
Karman’s spectrum. A random phase angle 7 is attributed to each Fourier 
mode. A final expression for the turbulent field is obtained as a sum over 
all the modes. 


Wave vector geometry of a single Fourier velocity mode 



- wavenumber k is randomly oriented: P(0) = sin 6/ 2 and P(<f>) — l/2rr 

statistical isotropy with respect to k 

- velocity vector a (k) -5> in plane perpendicular to k 

ip is random 

a(fc) ~ y^E(A:)Afc 

phase 7 of a (k) is random 


- «t(x) = Y I | COS (kj • x + 7) 
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A linear acoustic wave equation is derived. First and second order tur- 
bulence effects are retained in the wave equation. We assume that the wave 
propagation is lossless (isentropic) and that the turbulence is frozen. We 
only consider velocity fluctuations and do not include thermal fluctuations. 
e a is the acoustic Mach number and £t is the turbulence Mach number. 

A regular perturbation scheme in the turbulence Mach number is em- 
ployed to solve the wave equation. The N wave generated at the focus of 
the mirror is the boundary condition and it is represented as a sum of its 
Fourier components. 


2. Linear acoustic wave equation 

• lossless wave propagation 

• the turbulence is frozen 

• only turbulent velocity field is present 

• wave equation 

□=p=-4-ur-V|f-2^^-(;VpcJ() ( Mr 

^OO 

+ 2^4(/(ur'V)/V^ 

+ 2 l£ i 5i :(/ (/ Vpdt • V)u T dt\ - 4-(u t . V)(u r . Vp) ^r 2 

^OO 

-b-j— pV-((U7' • V)!^) 

C oo 

3. Perturbation scheme 

P = P0 + €tPi +€t 2 P2 + - 

B.C. p(0,t) = p N (t) 

Pn(0 = ^2 b n sin nu> 0 t 
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Results from the first order solution are shown here. Each graph repre- 
sents the sum of the incident N wave and the first order pressure perturba- 
tion after propagation through a single realization of the turbulent velocity 
field. As one observes, the first order pressure perturbation is responsible 
for the distortion of the N waves. Variations in waveform from peaked to 
rounded are noticed. Double-peaked and multiple-peaked waveforms are 
also shown. In some cases (e.g., the rounded waveform), the rise time of 
the waveform is changed, in other cases it is unaltered. However, the ar- 
rival time of each wave is the nominal arrival time (i.e., that of the incident 
wave). A calculation of the second order pressure perturbation is necessary. 
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The second order equation that has to be solved contains two secular 
terms and regular terms. A solution for the second order pressure pertur- 
bation is obtained. 

The total solution consists of the incident wave, the first order pressure 
perturbation, and the second order pressure perturbation due to the secular 
and regular terms. A renormalization technique is used to strain the z- 
coordinate. The straining of the z-coordinate is used to remove one of the 
singularities. The final solution is then written as a sum of the Fourier 
components of the N wave. The second order singularity introduces changes 
in the phase speed. At second order, the phase speed becomes dependent 
on the turbulence characteristics. 


5. Second order perturbation solution 

□ 2 £t 2 P 2 = (a - 0 - 7 )sinnw 0 (f - - 6cosnu 0 (t - , pa(0, t) = 0 

+ regular terms 

£t 2 P2 = — ~ — ^^zcosnw 0 (t- £-) + sin nu} 0 (t - 

7r)h>r> \ c ~/ 2 V 000 ' 


2 nu 0 

+ regular terms 


6. Total solution 


p = b n sin nu 0 (t - ^) + <-tP i 


(a -13- 7)c 0 


-zcos nu 0 




+ z sin nu 0 

2nu 0 


(<- i) 


+ regular terms of order e T 2 or higher 

Renormalization technique (strained coordinate) 
z = s(l + Ct 2 W 1 + • • •) 

Final solution P = 53 (* + sin nu 0 ^ t — 


• CtPi is first order perturbation solution 
2 (a - 0 - 7) 


• Ct Wi — — 








Results from the numerical simulations that incorporate the second or- 
der effects are shown. The phase speed is a function of the turbulence 
characteristics, and the actual phase speed is different from the nominal 
phase speed. Small variations in arrival time are observed. It is seen that a 
combination of first and second order effects of the turbulent velocity field 
is needed to fully explain the waveform distortion and the change in rise 
time. 


7. Results from numerical simulations 
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In order to make a convincing statement that the theoretical model is 
capable of simulating the propagation of plane N waves through a turbulent 
medium, we compared computed waveforms of the theoretical model with 
actual measured waveforms of the model experiment. The upper traces 
represent waveforms computed by the theoretical model, while the lower 
traces show waveforms from the model experiment. 

Two examples are given of a spiked waveform, and two examples of 
a rounded waveform. As one notices, the waveforms from the theoretical 
model exhibit the same distortion and change in rise time as that of the 
model experiment. 
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Comparisons of double-peaked waveforms axe shown, and also a com- 
parison of a U-shaped waveform and a rounded waveform is presented. 
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We showed that a model experiment can be successful in simulating the 
propagation of sonic booms through atmospheric turbulence. 

We also reviewed sonic boom models and compared the data from the 
model experiment with the results from the models. We found that only 
Pierce’s wavefront folding model is fairly accurate and that results from 
other models are not confirmed by the model experiment data. 

A new theoretical model is developed in which plane waves propagate 
through single realizations of a turbulent velocity field. The wave equation 
is solved by a perturbation method. The first order pressure perturbation 
creates the distortion of the N wave, and at second order a singularity 
occurs. The second order singularity introduces changes in the phase speed. 
The results from the theoretical model are confirmed by comparison with 
measured waveforms from the model experiment. 


Conclusion 

- model experiment is successful 

- reviewed sonic boom models 

only Pierce’s model is fairly accurate 
others are not confirmed by model exp. data 


- developed new theoretical model 

waves are propagated through single realizations of turbulent 
velocity field 

second order solution introduces a dispersion effect 
waveform distortion and change in rise time is caused by dispersion 
effect 

results are confirmed by comparison with measured waveforms 
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Slide 1 — 


This is the familiar formulation of sonic boom propagation. The near-field 
signature strength is defined by the F-function. There is an amplitude factor, 
the inverse of root-B, which is a generalization of cylindrical spreading. 
Nonlinear steepening appears as an adjustment to arrival time. This definition 
of the age parameter is very convenient, since once it's been computed for 
given flight conditions it can be applied to any F-function. 


SONIC BOOM AMPLITUDE AND AGING 
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Slide 2 — 


The kind of F-function we're most interested in is a shaped minimized one. 

This is George's F-function for a minimum-shock boom. (For simplicity. I'm 
only discussing the forward half, hence George's original form rather than 
the George/Seebass or Seebass/George extension to front and rear shocks.) As 

it ages, the initial delta-function impulse (Jones's asymptotic optimum) ages as 
an N-wave, while the isentropic compression behind follows. Everything is a 
very simple function of age parameter. 


GEORGE MINIMIZED F-FUNCTION 
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Slide 3 — 


The target ground boom occurs when the N-wave from the impulse just 
coalesces into the ramp, and the ramp still has some rise time which is slow 
enough to be not audible. Everything — including the matching value of the 
initial impulse — is related in simple ways. It’s worth looking at the age for 
which the ramp turns into a shock, as well as the age for the design condition. 
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These values are idealizations of the various low-boom designs which have 
been discussed over the past few years, for minimum boom Mach 1.6 to 1.8 and 
flight altitude of 40 to 50 kfeet. I've idealized the duration and slope of the 
isentropic compressions, and assumed a perfectly matched nose impulse. The 
design age parameter is about 0.8, with full shock coalescence occurring 
around 1.3. For this signature's initial and target ramp durations, if the 
signature ages by perhaps 30 to 60 percent more than the optimum age it will 
degenerate into a very noisy full shock. Lesser degrees of "overaging" will 
not increase its loudness as drastically. 

The 50 msec ramp is somewhat arbitrary. With Leatherwood, Sullivan, and 
Shepherd's excellent results from NASA — Langley's boom box experiments, it 
would be appropriate to establish formal target values for the ramp slope based 
on optimizing loudness. 


TYPICAL VALUES 

M = 1.7 

T 0 = 0.12 sec (L « 200 ft) 

= 0.06 Vft\ F 2 = 0.15 VF 
x = 0.05 sec 

Resultant Design Condition: 

A = 0.8 sec YfF" 1 

Shock Coalescence Condition: 

A shock = 13sec 
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Slide 5 — 


In a uniform atmosphere, the age parameter grows as the square root of 
distance. A convenient model for parametric analysis is an isothermal 
exponential atmosphere, with straight ray paths and fairly simple complete 
expressions. The real atmosphere can reasonably be approximated by a scale 
height in the range indicated. The isothermal-exponential atmosphere age 
parameter has an asymptotic limit, which equals the uniform atmosphere 
value at a radius (in the uniform atmosphere) which is about a scale height. 
This asymptote leads to the concept of freezing. 
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Slide 6 — 

Here are uniform and isothermal-exponential atmosphere age parameters at 
Mach 1.7. The asymptotic "frozen" limits are indicted on the right edge. Notice 

that these limits are not reached, or even approached, until very high 
altitudes — perhaps double the 40 to 60 kfeet of interest for HSCT. The term 
"freezing" has been used fairly often, and this familiarity has led to some 

common misconceptions. Atmospheric gradients clearly slow the aging 
process — to the point that McLean’s midfield signature concept is practical — 
but we are not close to the freezing regime. The isothermal-exponential age 

parameters in this figure, in the altitude range of interest, look like 
diminished versions of the uniform age parameter. 

Having the correct age parameter is essential even for predicting N-wave 

sonic booms. The overwhelming success of that endeavor suggests that there 
is no need to conduct elaborate experiments to prove that we know how sonic 
booms age, and certainly no need to become obsessed with asymptotic freezing. 

What certainly does need attention is determining whether mid-field 
signature aircraft designs are practical, and whether the signatures will 
survive under real-world conditions. Wind tunnel tests confirming 
configuration concepts have been successful, but wind tunnels are not large 
enough to allow aging to the midfield. Flight tests using RPVs and modified 
existing aircraft are being planned, and are necessary. The rest of this paper 
concentrates on proper scaling, with regard to age parameter, of these 
reduced-scale flight tests. Some phenomena which do not scale are identified, 

and emphasis is placed on pre-test analysis of elements which may not be 
intuitively obvious. 


Age Parameter — M = 1.7 



° 0 10 20 30 MO SO 60 70 80 

R , fc j 1 o 4 e e t 
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Slide 7 


This is an age parameter chart similar to the previous one, except that it is at 
Mach 2.0. Differences between this chart and the previous one exhibit the 
importance of using the correct Mach number. That the age parameter is 

bigger at Mach 2 than at 1.7 illustrates one of the reasons why boom 
minimization is easier at low Mach numbers. 

The age parameter for the real atmosphere is somewhat more complicated, 
depending on both the flight altitude and the ground altitude, rather than just 
the difference. Flight test design must use the real atmosphere. However, the 
isothermal-exponential model is adequate for the purposes of the current 
discussion. 


Age Parameter — M = 2.0 
Uniform and Isothermal Atmospheres 
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Slide 8 — 


A real-world consideration is that the age parameter is larger off track than 
on. This is the relation for the isothermal-exponential atmosphere. At an 
azimuth of 45 degrees, the age parameter is about 15 percent bigger, and at 
60 degrees it is about 40 percent bigger. There is a favorable benefit with the 
amplitude factor, but care must be taken to allow a margin before shock 
formation. In a flight test of an axisymmetric vehicle, on- and off-track 
measurements can be used to obtain several ages per flight. 


OFF-TRACK AGING 

• On-Track: <J> = 0, r = z. 

• Off-Track, Isothermal Model: 

- r = z/cos <(> 

- H effective "> H /cos $ 

A = V 2 71 H / cos f erf ( Vz/2H ) 

2 Y2p P a a 

\ = \=0 /V cos ? 
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Slide 9 — 


The first detail to consider for scaling flight test experiments is how the source 
scales with size. This is a simple representation of the volume and lift 
components of the F-function. The volume component always remains the 
same. Increasing the aircraft altitude at Fixed weight increases the 
F-function, as well as the age parameter due to increased distance. Flying at 
an optimum lift coefficient keeps F Fixed, and there is some recovery from the 
reduced amplitude factor at higher altitudes, so long as the increased aging 
does not lead to serious shock coalescence. This is why mission profiles with 
increasing altitude as fuel bums off tend to not show loudness increases, but 
rather decreases. 

For model scaling purposes, it is reasonable to assume that the boom is 
dominated by volume, or that full-scale lift coefficient will be replicated. The 

entire F-function then behaves like volume boom. 


F-FUNCTION FOR VOLUME AND LIFT (FIXED Ml 


F(x) 


j_ f* A" (|) ,, Pcos<|> i j-x L'(g) 

2lt J 0 (x- y>/2 ** P~uJ* J o (X- 5)1/2 

v v t V v , 


Volume 


Lift 


• Volume F-Function Fixed by Geometry. 

• Lift Component: 

- Increases if Raise Altitude at Fixed L 
(Off-Design Condition) 

- Decreases if L Decreases at Fixed Altitude 
(Fuel Bum-Off at Fixed Altitude) 

- Stays Constant at Fixed C L 
(Optimized Cruise Profile) 
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Slide 10 — 


For geometrically similar models, the F-function scales as the square root of 
length. A model less slender than full scale (larger f) increases F, which 
would in turn require a smaller propagation distance to retain the same 
relative aging. 


SCALING OF VOLUME BOOM 


Let 


F V M 2n 


r ^ 


A (x) = f 2 l 2 A (x/1) 


1 = length 
f = fineness ratio 
A (x/1) = nondimensional shape 
F v = f 2 fT F v (x/1) 
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Slide 1 1 — 


To achieve proper scaling in a uniform atmosphere simply requires altitude 
proportional to model length. In a real atmosphere, where the rate of growth 
of the age parameter is slower at larger distances, a model experiment requires 
flight altitudes which are smaller than obtained by using the model to full- 
scale ratio. If, for some reason, the model flight altitude must be proportional 
to model length, then the model must have a more slender fineness ratio. 


MODEL FLIGHT TEST 

• Require (AF) moc j e j — (AF)f u u scale 

• F oc yi/L , so need A y 1/L 

• Uniform Atmosphere: A °c i r- 
So Scale r model /r^,, - 1/L 

• Real Atmosphere: A tapers off 

^ r model/ r full < 
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Slide 12 — 


This table shows the scaling for three size models: full, one-third, and one- 

tenth. This roughly covers the range of proposed tests. The quantity shown as 
F is the nominal full-scale F-function at the end of the ramp, which was 0.15 in 
our example. The age parameter is that required so that age times F will scale 
proportionately to model length; for full scale, the age parameter for 50,000 
feet is shown. The value of z shown for the models is the corresponding 
altitude in the real atmosphere. The altitude ratio is always less than the model 
size ratio, and for the smallest scale is about half. The model age parameter is 
also proportionately more sensitive to altitude errors than full scale. It is 
straightforward to calculate required altitude precision, and the test plans 
must address this. 

The final parameter in the table is the duration of the ramp. A nominal 
50 msec full-scale ramp is expected to be clearly discernible. The corre- 
sponding model ramp durations are very close to the 3 to 10 msec range 
normally encountered for shock wave rise times. If we expect shocks to have 
their typical "full-scale N-wave" rise times, the smaller scale tests are 
somewhat dubious. 


RELATIVE 

SCALES 

FOR ONE-THIRD 

AND 


ONE-TENTH 

MODEL 

TESTS 


Scale 

F 

A 

Z 

z/z fiIll 

T 

1.0 

0.15 

0.84 50,000 

1.0 

50 msec 

1/3 

0.087 

0.48 

9,500 

0.19 

17 msec 

1/10 

0.047 

0.27 

2,600 

0.05 

5 msec 
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Slide 13 


Rise times may not, however, be as long as they are at full scale. Bass, Layton, 
and Bolen ( JASA , July 1987) measured projectile shock waves which were 
considerably thinner than expected from steady theory. The reason was that 
propagation distances were sufficiently short that they had not had time to 
reach steady state. At last year's sonic boom workshop, Raspet discussed 
"healing times" and suggested that even full-scale N-wave shocks might not be 
steady. If full-scale shocks have not achieved steady thickness, then model 
shocks must certainly be investigated. 

The current HSR team has developed several models which are suitable for 
calculating the evolution of shock structures. This is the time to use them. 


POTENTIAL FOR THIN SHOCKS 

• Unsteady Shock Formation Over Short Distances 
(Bass, Layton, and Bolen - JASA, July 1987). 

• Healing Time of Weak Shocks (Raspet, 1992). 

• Current Analysis Tools Available to Predict. 
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Slide 14 — 


As a final element, I'd like to show some "anything can happen" data. This is 
of interest because this could mask results of an "aging" flight test, and also 
because it provides an example of turbulence effects which do not fit into the 
usual expected pattern. 

We recently performed what should have been a very dull Boomfile test 
protocol on an aircraft not yet in the data base. Tests were conducted in the 
early morning, under overcast skies, so there was no convective turbulence. 
There was virtually no surface wind. The weather had been changing 
(it snowed overnight, and was clearing during the tests) and a rawindsonde at 
a site about 100 miles away indicated wind shears at about 15 to 20 kft and 20 to 
25 kft MSL. The site was at about 5 kft MSL. 


SUPPLEMENTAL BOOMFILE FLIGHT TEST 

• Documentation of Steady Boom for an Aircraft 
Not Previously Measured. 

• Flight Parameters: 10k, 18k, 25k AGL, 

at Mach 1.2, 1.25, 1.30. 

• Recorded On and Near Centerline by USAF 
BEAR Systems. 

• Early Morning, Calm Surface Conditions. 

• Rawindsonde During Test, About 100 Miles Away. 
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Slide 15 — 


Here is a typical good measurement. In this test, which yielded about 25 valid 
recordings from ten passes, about half the records were of this quality. The 

rest had significant distortion. The next three slides are particularly 
interesting. 


25 kft, M = 1.3 


File wO‘l 1930H. 130 00:00:00.00 

Pi»ax= M9 Pmin = -.96 1801 points Site 0 S/N 0 

105.5 d8C 1 15. 7 dB (flat) 89.2 dflfl 103.9 PLdB 
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Slides 16, 17, 18 — 


These are distorted measurements at 25, 18, and 10 kfeet AGL, respectively. 
Overpressures and durations were as expected. There was considerable 
distortion, and rise times were very long. These three booms have a particular 
characteristic that they look very much like minimum-shock shaped booms. 
Slides 15 and 16 are the same flight condition, so the nice minimum shock 
shape in 16 is anomalous. 

Through most of this paper I have discussed simulation issues that could lead to 
what the statisticians call Type I errors — failing to confirm a phenomenon 
that exists. These three booms raise the possibility of Type II errors — 
detecting an effect that is not there. 

These measurements are also interesting because they have features not fully 
consistent with the usual characteristics of distorted booms. Distortion was 
more common for the higher altitude runs, while the low-altitude runs yielded 
a higher percentage of clean booms. Surface layer turbulence would affect all 
altitudes. There was a considerable amount of noise after the booms — the kind 
that sounds like the echoing of distant thunder. This tended to persist for 
several seconds. Each BEAR would typically record three or four records per 
boom, with the first being the N-wave and the rest being noise. This was much 
more than observed during the 1987 Boomfile tests at Edwards AFB. 

The site was flat, so the distortion and aftershocks had to be a combination of 
scattering and multiple paths. The long aftershocks and the substantial 
distortion are consistent with scattering from the shear layers and weather at 
higher altitudes, i.e., not in the surface and mixed layer as normally expected. 
Anecdotal reports of "echoing" in NASA's JAPE II sonic boom propagation 
experiment should be re-examined in terms of multi-path propagation from 
higher altitude atmospheric structures. 

These measurements underscore that individual test conditions can obscure 
fine details such as would be seen in model tests of midfield signatures. Both 
Type I and Type II errors are possible. Care must be taken to design the initial 
round of shaped boom flight tests so as to avoid these conditions. It is also clear 
that atmospheric conditions occur which lead to distortions somewhat 
different from those seen in summertime desert tests. It would be prudent to 
extend propagation investigations into other geographical areas and 
atmospheric conditions. 
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12.5 kft, M = 1.2 


File utO'1201 7fl. 130 20:17:15.18 January 30 1993 

P m ax= 2.21 Pmin = -1.80 1266 points Si te '1 S/N >1019 

106.2 dBC 121. 1 dB I flail 88. 1 dBR 102. ‘I PldB 
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The purpose of this paper is to summarize the latest three sonic boom 
laboratory studies performed at NASA-Langley using the Sonic Boom Simulator. 
The first used synthesized idealized outdoor boom shapes which were filtered to 
represent booms heard inside a house. The test explored the efficacy of various 
metrics in assessing both loudness and annoyance responses to these booms. 
The second test investigated the effects of adding single reflections to idealized 
boom signatures, and the third compared booms recorded from real aircraft with 
idealized boom signatures to determine if subjects rated the real booms differently. 
In these studies, as in previous studies performed at NASA-Langley, there was a 
continuing effort to evaluate metrics for predicting the subjective effects of sonic 
booms. 


OBJECTIVES 


• Summarize recent lab studies 

— Quantify indoor/outdoor subjective effects 
— Determine effects of ground reflections 
— Quantify loudness response to recorded booms 
— Evaluate metrics 
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The tests were conducted using the NASA-Langley sonic boom simulator. 
The simulator is a person-rated, airtight, loudspeaker-driven booth capable of 
accurately reproducing user-specified waveforms with peak pressures of up to 
138 dB. Signals are preprocessed to compensate for non-uniformities in the 
frequency response of the booth and sound reproduction system. The system 
is fully described in reference 1. Although tests using the Simulator cannot 
completely replicate conditions in real life, they allow listeners to evaluate sonic 
booms under controlled conditions. 


SONIC BOOM SIMULATOR 
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A total of 168 subjects took part in the three tests: 72 in test 1, 48 in test 2 
and 48 in test 3. Tests 2 and 3 were run concurrently and used the same subjects. 
Magnitude Estimation Scaling was the psychometric methodology used. In this 
methodology, one boom was selected as a reference and given a score of 100. 
Subjects were asked to compare all other booms to this reference on a ratio basis. 
The wave shapes used in the tests included simulations of the N-waves, front- 
shock minimized (shaped) booms, composite booms (that is, an original N-wave 
or minimized boom waveform, to which was added a single delayed copy of the 
original), and recorded booms, which were digitally recorded and modified for 
reproduction in the simulator. 


EXPERIMENTAL APPROACH 


• Three experiments 

• Magnitude estimation scaling 

• 168 test subjects 

• Signature shapes 

— N-waves 
— Shaped 

— Composite (original + reflected) 
— Recorded 
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Five metrics have been studied in these and previous tests run at Langley. 
These were unweighted sound exposure level (SEL), A-weighted sound exposure 
level (SEL a ), C-weighted sound exposure level (SELc), perceived level (PL) 
using the Stevens’ Mark VII procedure, and Zwicker loudness level. Considering 
the results of all tests run in the Simulator to date, PL has proved to be the most 
effective metric for predicting subjective reactions to sonic booms. Therefore the 
results of the three tests in the present paper are reported in terms of PL. Full 
details on the calculation procedures are given in reference 2. Previous tests are 
described in references 3, 4 and 5. 


METRICS 


• Five metrics considered 

— SEL (unweighted) 

— SELA 
— SELC 

— Perceived Level (Stevens’ Mark VII) 

— Zwicker Loudness Level 

• All studies have evaluated metric performances 

• Perceived Level selected as best metric for general use 
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The objectives of the first experiment were to investigate whether people 
respond differently to sonic booms heard indoors as compared with booms heard 
outdoors, to determine whether there would be a difference if they were asked 
to rate loudness or annoyance, and to validate the PL metric. Tests in the booth 
cannot be true simulations of reactions in the home, as there is no rattle, no fear of 
property damage, and the subjects’ sole activity is to sit and listen to the booms. 
Hence this test could only compare loudness and annoyance responses in the lab. 
The test results do provide a basis for determining which descriptor was more 
appropriate and for evaluating which metrics predicted people’s reactions most 
accurately. This test used two groups of 36 subjects each; one group was asked 
to rate the booms based on annoyance and the other was asked to judge loudness. 
Both groups heard the same signals played at the same levels; both used the same 
reference boom with a score of 100. The calculated metric levels were unchanged 
between the two parts of the test; only the subjective ratings differed. 


EXPERIMENT 1 OBJECTIVES 

• Quantify loudness and annoyance of simulated indoor and 
outdoor booms 

• Determine appropriate subjective criterion measure 

• Verify PL metric 
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The signatures used in this study were based on three shapes: an N-wave, 
a front-shock minimized boom with a ratio of front shock to peak overpressure 
of 0.75 and a front-shock minimized boom with a ratio of front shock to peak 
overpressure of 0.5, as shown in the figure. The front-shock minimized (FSM) 
booms had a secondary rise time of 60 msec. All booms had a duration of 300 
msec, and for all three shapes front-shock rise times of 2, 4 and 8 msec were 
used. Thus a total of nine waveforms comprised the outdoor signals. 


NOMINAL OUTDOOR SIGNATURES 


(All durations = 300 msec) 
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Each outdoor shape was digitally filtered using two filters, one representing 
the attenuation through a typical house wall with windows open and the other with 
windows closed. Both gave no attenuation below 10 Hz, which is considered the 
simplest reasonable assumption in view of the lack of data at these frequencies. 
Above 60 Hz, these shapes were mathematical curves based on measured data. 
For both filters, the effect was to remove high frequency energy, more of which 
was removed by the “windows closed” filter than by the “windows open” filter. 


HOUSE FILTERS 
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This chart shows simulated indoor and outdoor booms, as measured within 
the simulator, for cases where the outdoor booms had rise times of 2 msec. 
Differences between the N-wave and the FSM booms can still be seen in the 
filtered shapes, but the high frequency, sharp corners of the originals have been 
greatly reduced. The nine original shapes were each presented to the subjects in 
the original form (the “outdoor” booms) as well as in the two “indoor” forms 
(“windows open” and “windows closed”). Measurements of sonic booms inside 
real houses would show more complex signatures, because of phase changes which 
were not simulated. However, these changes do not alter the frequency content 
of the boom. No reverberation effects were included in these simulations. 


OUTDOOR AND INDOOR SIGNATURES 

(rise time = 2 msec) 





N-WAVE 


MIN50 


MIN75 


1G1 



LOG{GEOMTrRIC MEAN) 


Results obtained in the first experiment are summarized in this chart which 
displays plots of the logarithm of the geometric means of the subjective ratings 
versus Perceived Level. The geometric mean is the appropriate central tendency 
measure for data obtained using the Magnitude Estimation psychometric method. 
For the outdoor booms, the results (displayed in the plot on the left) show no 
difference between loudness and annoyance responses, which agrees with previous 
studies. However, differences between loudness and annoyance scores were found 
for the indoor booms. For both simulated indoor conditions, the annoyance 
responses where higher than loudness responses. The plot on the right shows 
the results for the “windows closed” condition. The “windows open” results were 
similar, except that the differences between loudness and annoyance responses 
were smaller. This is reasonable since the “window open” filter had less effect 
on the boom than the “windows closed” filter. Statistically the slopes are not 
significantly different for either pair of lines. At the levels used in this test, the 
low frequencies were dominant in the indoor booms. Thus it could be inferred that 
they caused an increase in annoyance greater than their contribution to loudness. 
Based upon these results it is recommended that studies of indoor booms use 
annoyance rather than loudness as the judgement criterion, as the use of loudness 
may underestimate the booms’ unacceptibility. 


LOUDNESS AND ANNOYANCE RESPONSE 


OUTDOOR 



INDOOR-WINDOWS CLOSED 
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The performance of PL for judgements of loudness (left plot) and annoyance 
(right plot) for each simulated listening condition is shown in this figure. For a 
given PL value, the indoor booms were rated lower in loudness than the outdoor 
booms. The “windows open” condition lies between the outdoor and the “windows 
closed” loudness results. However, the indoor and outdoor booms of the same PL 
level were rated equally annoying. In both plots, the slopes of the regression lines 
for the indoor boom results differ significantly from those for the outdoor boom 
results. Results for SELA or the other metrics investigated (not shown) clearly 
indicate that PL predicted more accurately the annoyance for the combined data 
set of indoor and outdoor booms. 


PL METRIC PERFORMANCE 



ANNOYANCE 



WINDOWS CLOSED 


72 


77 82 67 92 

PERCEIVED LEVEL, dB 


163 



Field measurements of sonic booms are usually made with flush-mounted 
microphones, in which case a single reflection is present in the recording with 
a zero time delay between the direct boom and the reflected boom. However, 
outdoor listeners will usually receive a direct boom followed by a reflection off the 
ground having a finite delay, of the order of 8 msecs. In the second experiment, 
the objective was to evaluate the effects of different delay times on subjective 
response and to find a metric that accounted for these effects. To investigate this, 
single “reflections” were added to idealized boom shapes. Delay time between 
the “direct” boom and the “reflected” boom was a variable in this study. As only 
outdoor boom shapes were used, subjects were asked to rate loudness. 


EXPERIMENT 2 OBJECTIVES 


• Quantify loudness of sonic booms containing reflections 
(direct + single reflection) 

• Determine delay time effects 

• Verify PL metric 


This chart demonstrates the effect of adding to an idealized N-wave shape (the 
“direct” boom) a single time-delayed version of the same shape (the “reflected” 
boom). The combined wave form is called a “composite” boom. No phase 
change or attenuation was introduced into the delayed wave. If the delay is 
zero, the result of combining the direct and reflected booms is a doubling of 
overpressure while the rise time remains unchanged. If the delay is equal to the 
rise time, the resulting wave has a rise time twice that of the original, while the 
overpressure is nearly double the original. Other delay times result in composites 
having complex, multi-segmented pressure increases. 


COMPOSITE BOOM SIMULATION 


DELAY TIME 



NO DELAY DELAY - RISE TIME DELAY > RISE TIME 





DIRECT + REFLECTED = COMPOSITE 
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The basic shapes considered were an N-wave and a front-shock minimized 
boom with a ratio of front shock overpressure to peak over pressure of 0.5. 
For both shapes rise times of 3, 6 and 9 msec were used, which resulted in six 
waveforms for the “direct” booms. For each direct boom, six values were selected 
for the delay time. Delays of 0 msecs and 12 msecs were used for each direct 
boom. A value of delay equal to the front shock rise time was also selected, 
together with rise time + 1 msec and rise time — 1 msec. The sixth value of 
the delay was 3 msecs, except for the booms with 3 msec rise times, where a 
value of 8 msec was used. 


EXPERIMENTAL WAVEFORMS 


• Six direct boom waveforms 

— 3 N-waves (rise times = 3, 6, 9 msecs) 

— 3 shaped (front shock rise times = 3, 6, 9 msecs) 

• Six reflection delay times for each direct boom 
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The composite waveforms that resulted when a 6-msec N-wave was used 
as the original waveform are illustrated here. In these idealized plots, the direct 
boom has an overpressure of 1 psf. The delay time between the “direct” boom 
and the “reflection” was varied between 0 and 12 msec. The zero delay results in 
a waveform with a rise time of 6 msec, but a peak overpressure of 2 psf. When 
the delay equals 6 msec, the composite wave has a peak overpressure of almost 
2 psf, but the rise time has been doubled to 12 msec. The other delays produce 
other more varied waveshapes. 


COMPOSITE WAVEFORMS 


(N-wave, T = 6 msec) 
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The effects of varying the delay for booms with a range of rise times were 
predicted using the PL metric. This chart shows the predicted PL values for N- 
waves with rise times ranging from 0.1 msec to 12 msec, and delays from 0 to 14 
msec. The curves show a maximum at zero delay, and a minimum for most cases 
when the delay equals the rise time. Some of the irregularities in these curves 
are due to the fact that they were calculated using 1/3 octave bands. As the delay 
time changes, there are subtle changes in the spectrum, which cause energy to 
shift in frequency. As it shifts from one band to another, there are changes in the 
calculated values. If a continuous algorithm is used, such as that used for SELA, 
these irregularities are removed, SELA shows the minimum when the delay 
exactly equals the rise time, but for the short rise times PL shows the minimum at 
somewhat greater delays. This is probably an artifact of the calculation procedure. 


PREDICTED DELAY TIME EFFECTS 
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This chart shows the results obtained for the six basic waveforms, each with 
six different delay times. The solid lines show the subjective responses, converted 
to equivalent PL, while the dashed lines show PL calculated from the measured 
data. The results are normalized to 0 dB at a delay time of zero. In general there 
is good agreement between PL calculated from measurement and PL derived from 
the ratings. The ratings all show a minimum when the delay time is equal to the 
boom front shock rise time. Predicted PL and PL based upon measurements also 
gave minima when delay time equals front shock rise time except for the 3 msec 
booms which have a minimum at a 4 msec delay. 

Most sonic boom measurements are made with flush-mounted microphones, 
for which there is zero delay between the direct wave and the ground-reflected 
wave. Hence, these measurements will yield loudness levels that are conservative 
compared to outdoor situations, where the listener will experience a ground 
reflection with a delay typically of the order of 8 msec. In an enclosed situation, 
there could be more than one reflection, and more complex wave forms could 
result. 


DELAY TIME EFFECT 


TAU - 1 ma 


TAU ■ 6 ma 


TAU - 1 ma 





N-wavas 




Shapad Booma 
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LOG(GEOMETRIC MEAN) 


The figure shows a plot of subjective ratings as a function of PL calculated 
from measured data for the set of composite booms containing reflections with 
zero delay and for the set containing reflections with non-zero delays. There is 
no statistically significant difference between the two lines, and PL can be said 
to account for the effects of delay on loudness. 


DELAY EFFECT - PL METRIC 



PERCEIVED LEVEL, dB 


170 


The third experiment used real sonic boom recordings made at White Sands 
Missile Range from T38 and FI 5 flyovers. The objective of the test was to 
determine if people reacted differently to real booms as compared to idealized 
booms, and to assess metric performance for predicting the loudness of real booms 
that have been distorted by propagation through the atmosphere. The recordings 
were made using the United States Air Force’s BEAR (Boom Event Analysis 
Recorder) systems using a 8 kHz sample rate. Seven hundred recorded signatures 
were scanned and thirteen booms selected that fitted into four major categories 
of sonic boom shape: 

1. N-wave 

2. Rounded 

3. Peaked 

4. U-shaped 

The recorded booms were digitally edited to remove some of the background noise 
and to increase the time between their front and rear shocks to 300 msec, which is 
more representative of a HSCT. The front and rear shocks were not altered by this 
procedure. Three idealized N-waves and two boom shapes based on predictions 
for HSCT from various CFD codes were also included, resulting in a total of 
18 signals. The waveforms based on the CFD predictions were simplified, but 
both contained more than two pressure peaks. Because of the existence of extra 
shocks between the front and rear shocks, these were designated “intermediate” 
booms. All the signals, recorded and idealized, were preprocessed to account for 
the booth frequency response. 


EXPERIMENT 3 OBJECTIVES 


• Quantify loudness of recorded signatures 

• Compare with results for idealized booms 

• Verify PL metric 
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This chart presents some of the results from this test, together with three 
examples of the recorded waveforms as reproduced in the simulator. These 
data show that responses to the four categories of real booms did not differ 
from responses to the idealized N-waves when plotted against PL calculated from 
measured signals. Thus PL accounted for any differences between the waveforms. 
These booms were not played at the same levels as they were heard at White 
Sands, nor even at the same relative levels; instead they were adjusted to cover 
the same range of SELa values. 


RESPONSE TO WHITE SANDS BOOMS 
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The only category of boom that stood out in this study were the “intermediate” 
booms. Time histories of these waveforms measured in the simulator are shown on 
the right of this chart. The plot displays the same judgement data as the previous 
chart, but the results for the two intermediate booms have been added separately. 
The intermediate boom results fall somewhat below those for the other boom 
categories, indicating that the subjects rated them as having lower loudness than 
the PL metric predicted. The slopes of the regression lines for the intermediate 
boom results differ from those for the other booms. Boom asymmetry is a possible 
contributing cause, though a previous study on asymmetry (reference 3) would 
predict little effect for Intermediate 1 and none at all for Intermediate 2. The 
multiple peaks may have some masking effect on each other. Further study is 
needed to understand these results, which were based on only two samples. 


RESPONSE TO INTERMEDIATE BOOMS 



PERCEIVED LEVEL, dB 



INTERMEDIATE 1 



INTERMEDIATE 2 
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In summary the findings of the three studies showed that loudness and 
annoyance gave equivalent results for outdoor booms but differed for indoor 
booms. Annoyance is therefore recommended as the most appropriate criterion for 
indoor boom judgements. For booms containing a single reflection, the loudness 
ratings were higher when the delay between the direct and reflected booms was 
zero than when the delay was greater than zero . Hence loudness values calculated 
from flush-mounted microphone measurements will be conservative. Subjects did 
not judge recordings of real booms any differently from idealized booms. The PL 
metric was validated for annoyance ratings, for indoor and outdoor booms, for 
booms with reflections and for real booms compared with idealized simulations. 


SUMMARY 

• Loudness and annoyance are: 

— equivalent for outdoor booms 

— not equivalent for indoor booms 

• Annoyance most appropriate criterion 

• Booms with single reflections are 

— less loud for nonzero delay times 

— Hence measurements with flush-mounted microphones 
are conservative 

• Real booms not judged differently from idealized booms 

• PL metric validated 
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Comparison of Methods of Predicting Community Response to 
Impulsive and Nonimpulsive Noise 


Abstract 


Several scientific, regulatory and policy-coordinating bodies have developed 
methods for predicting community response to sonic booms. The best known of 
these is the dosage-response relationship of Working Group 84 of the National 
Academy of Science’s Committee on Hearing, Bioacoustics and Biomechanics 
(Galloway, 1981). This dosage-response relationship between C-weighted Day- 
Night Average Sound Level and the prevalence of annoyance with high energy 
impulsive sounds was derived from limited amounts of information about 
community response to regular, prolonged, and expected exposure to artillery and 
sonic booms. 

U.S. Army Regulation 201 adapts this approach to predictions of the acceptability 
of impulsive noise exposure in communities. This regulation infers equivalent 
degrees of effect with respect to a well known dosage-response relationship for 
general (nonimpulsive) transportation noise. Differences in prevalence of 
annoyance predicted by various relationships lead to different predictions of the 
compatibility of land uses with sonic boom exposure. An examination of these 
differences makes apparent several unresolved issues in current practice for 
predicting and interpreting the prevalence of annoyance due to sonic boom 
exposure. 
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BACKGROUND 


Any systematic approach to predicting and interpreting community response to 
noise exposure requires solutions to four fundamental problems: 

1) Definition of community response; 

2) Characterization of noise exposure; 

3) Derivation of a predictive relationship between "community response" 
and noise exposure; and 

4) Inference of regulatory policy from one or more predictive relationships 

Finding practical answers to these questions requires detailed attention to issues 
such as the ease and cost of measuring selected quantities, the desired accuracy 
and precision of predictions, and relationships among alternate metrics of noise 
exposure, community response, and land use compatibility. 

Most of the debate about compromises and assumptions needed to predict and 
interpret community response to transportation and other non-impulsive noise was 
conducted in the 1970s. To make a long story short, a collection of federal 
agencies (FICUN) adopted a common approach - based in large part on the work 
of Schultz (1978) - built on the following assumptions: 

1) that "community response" can be usefully treated for most purposes 
(and in particular, for the airport neighborhood case) as the proportion of 
a residential community annoyed to a consequential degree by noise 
exposure 1 ; 

2) that a cumulative measure of outdoor A-weighted sound levels in 
residential neighborhoods incorporating a so-called nighttime penalty (the 
Day-Night Average Sound Level, or DNL, as developed in the 
Environmental Protection Agency’s 1974 "Levels Document"), suffices for 
characterizing noise exposure; 

3) that a dosage-response relationship between the prevalence of 
annoyance and DNL, derived from a curve fitting exercise, is adequate for 
predictive purposes; and 

4) that regulatory policy can be based on interpretations of land use 
compatibility made in acoustic terms alone. 

The basic notion in any dosage-response relationship is that whatever quantity is 
plotted on the abscissa as the predictor variable is uniquely responsible for 
whatever quantity is plotted on the ordinate. In this case, the proportion of a 
community highly annoyed by noise is assumed to be determined not by 
individually notable noise events, but solely by some integration of outdoor 
neighborhood noise levels over a prolonged period of time. 
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"Land use compatibility guidelines” were subsequently developed by interpreting 
predictions about the prevalence of annoyance derived from the dosage-response 
curve. "Compatibility" was treated as an issue of noise exposure rather than one 
of noise effects, as though land use compatibility were somehow a property of 
noise exposure per se. 2 

The compromises and assumptions of the 1970s are not the only ones that could 
have been made, nor are they necessarily the most appropriate for all purposes. 
However, they provide the basis for the most widely understood and applied 
approach to assessing impacts of non-impulsive noise exposure on communities. 

DERIVATION OF DOSAGE-RESPONSE RELATIONSHIPS 

The data set on which Schultz based his original (1978) synthesis of the community 
noise literature contained 161 data points. One recent update of the so-called 
"Schultz Curve” (Fidell, Barber and Schultz, 1991) is based on more than 400 data 
points, while FICON (1992) has developed another dosage-response relationship 
based on a subset of these points. Each data point represents a field observation 
of a pairing of an exposure value and a percentage of social survey respondents 
describing themselves as highly annoyed. The percentage of social survey 
respondents is treated as an index of the stable, steady state prevalence of a 
consequential degree of long term annoyance in the community at large. 

The utility of Schultz’s approach to assessing community response to non-impulsive 
noise exposure quickly led to efforts to apply similar methods to the case of high 
energy impulsive noises 3 . Working Group 84 of the Committee on Hearing, 
Bioacoustics and Biomechanics (Galloway, 1981) of the National Research Council 
of the National Academy of Science made the initial (and still best known) effort 
to adapt the methods developed by Schultz (1978) to the case of impulsive noise. 

CHABA Working Group 84 preserved Schultz’s definition of community response; 
modified reliance on DNL as a noise metric only to the extent of substituting C- 
weighted for A-weighted sound levels 4 ; and developed a dosage-response 
relationship from an eyeball fit to a small number of social survey observations 
about the annoyance of impulse noise. 5 

The resulting dosage-response relationship for impulsive noise (Galloway, 1981), 
illustrated in Figure 1, is as follows: 

% Highly Annoyed = 100/(1 + exp(11.17- 0.153LQJ,,)) Eq. 1 

where Lq^ is the C-weighted Day-Night Average Sound Level created by sonic 
booms. 

The form of this dosage-response relationship is a sigmoid given by a logistic fitting 
function. The sigmoidal shape is a reasonable one, given the need for asymptotes 
in the relationship in the vicinities of 0 and 100%. The prediction equation 
reflects a negotiated consensus of engineering judgments, and is intended as an 
approximate curve fit rather than the product of a formal statistical analysis. 
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INTERPRETATIONS OF DOSAGE-RESPONSE RELATIONSHIPS 


The CHABA Working Group 84 relationship is derived from a considerably 
smaller data set than any of the relationships for non-impulsive noise. Of a total 
of fourteen data points, five represent observations of the annoyance associated 
with exposure to artillery fire, while the remaining nine are all derived from the 
only extensive study ever conducted of community response in an urban area 
subjected to sonic booms over a prolonged period. 

For six months in 1964, residents of this city were exposed to a maximum of 8 
booms a day at overpressures of 1 to 2 psf. Figures 2 and 3 are linear regressions 
between the prevalence of annoyance and A-weighted and C-weighted exposure 
values, respectively, for this data set, as described by Galloway (1981). The former 
(C-weighted) regression accounts for 94% of the variance in the annoyance data, 
while the latter (A-weighted) regression accounts of 87% of the variance. 

The circumstances of impulsive noise exposure in the artillery and Oklahoma City 
studies which are summarized in the CHABA relationship are noteworthy: they 
were all familiar, expected, predictable and of long duration. In the case of 
artillery noise, respondents were residents of neighborhoods near fixed firing 
points. Daily artillery noise was a familiar part of the noise exposure environment 
in these respondents’ neighborhoods. In the case of the sonic boom exposure in 
Oklahoma City, advance schedules for the numbers and times of occurrence of 
sonic booms were well advertised, and the aircraft producing the booms flew a 
single flight track for the entire half year study period. 
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The most common form of community exposure to sonic booms today is that 
produced in the vicinity of military supersonic operations areas. Caution is 
required in predicting community response and land use compatibility in such 
circumstances. Sonic booms are generally experienced at unpredictable and 
relatively infrequent times as short duration daytime noise intrusions of widely 
varying level. There is no agreement comparable to that for the case of urban 
noise about the most useful approach to predicting the annoyance of this type of 
noise exposure. 

Use of DNL or CDNL to predict the prevalence of annoyance is based on the 
"equal energy hypothesis". The equal energy hypothesis expresses the notion that 
the number, level and duration of noise events are fully interchangeable 
determinants of annoyance as long as their product (energy summation) remains 
constant. In other words, quantification of noise exposure in DNL for purposes of 
predicting annoyance reflects a tacit theory: that people are indifferent between 
the annoyance of small numbers of very high level noise events of short duration 
and the annoyance of large numbers of compensatingly lower level noise and/or 
longer duration noise events. 

This hypothesis is the underpinning of a convenient method for measuring noise 
exposure for purposes of predicting annoyance. When used as a predictor variable 
in a dosage-response relationship such as that synthesized by Schultz, DNL 
accounts for about half of the variance in a set of field observations about the 
annoyance of general transportation noise. This demonstrates that the equal 



Figure 2: Linear regression between C-weighted sonic boom exposure and 

prevalence of annoyance in Oklahoma City (data from Galloway, 1981 ) 
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energy hypothesis can provide a useful account of the data over a range of at least 
20 dB, from values of about 55 to 75 dB. 


There is, however, little empirical evidence on which to base extrapolations of 
predictions of annoyance at the low values of CDNL associated with infrequent 
exposure to sonic booms. For example, all of the data summarized by Galloway 
(1981) at low values of CDNL represent reactions to artillery fire, not sonic 
booms. 

ASSESSING LAND USE COMPATIBILITY WITH IMPULSIVE NOISE 
EXPOSURE 

The notion that environmental noise impacts can be construed for global purposes 
in terms of land use compatibility may be traced through a chain of noise metrics 
and prediction methods four decades long to the pioneering work of Rosenblith 
and Stevens (1953). The latest embodiment of this approach may be found in the 
Appendix to ANSI Standard S12.40-1990. 

The U.S. Army has adopted the clearest guidelines among federal agencies for 
land use compatibility with high energy impulsive noise exposure. Chapter 7 of 
U.S. Army Regulation 200-1, "Environmental Noise Abatement Program" (dated 
23 April 1990) addresses land use compatibility issues with respect to impulse noise 
exposure. The recommendations in this regulation are based on an equivalence of 



Figure 3: Linear regression between A-weighted sonic boom exposure and 

prevalence of annoyance in Oklahoma City (data from Galloway, 1981) 


183 



exposure of C-weighted to A-weighted cumulative exposure units. This 
equivalence is made not in terms of exposure levels, but rather in terms of the 
prevalence of annoyance implied by the FICUN (1980) and ANSI S12.40-1990 
guidelines. 

Equivalent exposure is inferred in the Army regulation through comparisons of 
annoyance predictions of Schultz’s original (1978) dosage-response relationship for 
general transportation noise and CHABA’s (Galloway, 1981) dosage-response 
relationship for high energy impulsive noise. Table 1 compares the prevalence of 
annoyance predicted by Schultz’s and other relationships. 

Table 2 examines the implications of drawing land use compatibility inferences on 
the basis of these equivalences. In the range of interest, differences in the slopes 
of the logistic fitting functions described by Galloway (1981) give rise to an 
approximate 5 dB difference in the criterion levels as expressed in A- and C- 
weighted cumulative exposure units. The net effect is to increase the difference 
between A- and C-weighted units associated with the same prevalence of 
annoyance. Thus, whereas Army Regulation 200-1 relies on the dosage-response 
relationship developed by Schultz (1978) to equate an A-weighted DNL of 65 dB 
with a C-weighted DNL of 62 dB, the equivalence developed from FICON’s 
(1992) dosage-response relationship is to a C-weighted DNL of 60 dB. Even 
greater differences are apparent if the analysis is restricted to sonic boom data 
only. 


Table I. Percentage of Community Highly Annoyed ("%HA") Predicted by Several 
Dosage-Response Relationships over a Range of DNL/CDNL Values. 


FITTING FUNCTION 

sill. 
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Hunn 
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%HA ~ 0.8553L d „ - 0.0401L dn * -t- 0.00047L*, 5 
(Schultz (1978)1 

3.9% 

8.5 

153 

24.6 

36.9 

% HA " 100/(1 + «p(10.43 - .132L dn )} 
[USAF logistic fit to Schultz data] 

4.0% 

IS 

13.6 

23.3 

37.0 

% HA - 0.0360L dn l - 3.26451.,*, + 78.92 
[Fidelt, Schultz, and Barber (1991)J 

8.3% 

12.7 

18.8 

26.8 

36.6 

% HA - 100/(1 + otp(11.13 - .HIL^)) 

| [USAF logistic Qt to subset of Fidell, Schultz and Barber (1991)] 

3.3% 

6,5 

12 3 

22.1 

36.5 

| %HA - 100/(1 + ccp(11.17 - .1531^)) 
fj (CHABA WO 84 [Galloway, 1981]) 

2.9% 

a 

22.7 

38.7 
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Table 3 shows some examples of predicted consequences of several sorts of sonic 
boom exposures. Note that even small numbers of relatively modest booms (for 
example, four booms at 1 psf per day) can lead to a prediction of noise exposure 
inconsistent with single family residential use. 

Although there is little alternative to relying on such predictions on an interim 
basis, a larger body of direct evidence about the annoyance of sonic boom 
exposure is clearly needed. It is possible, for example, that the time constants of 
arousal and decay of annoyance with impulsive noise exposure may differ from 
those for non-impulsive noise. People living near future overland supersonic flight 
corridors might react more quickly or more vigorously to the novel experience of 
sonic boom exposure than to more familiar forms of noise exposure. 

CONCLUSIONS 

Confidence in predictions and interpretations of the effects on communities of 
high energy impulses in general and sonic booms in particular is not as great as for 
the case of general transportation noise. Small differences in assumptions and 
procedures may lead to large differences in assessments of the effects of impulsive 
noises on exposed populations. Some of these assumptions and procedures that 
have not been revisited for a decade could benefit from further scrutiny. 

The circumstances of noise exposure produced by sonic booms of en route aircraft 


Table II. Land use compatibility guidance inferred from equivalent prevalence of 
annoyance for A-weighted and C-weighted Day-Night Sound Levels for alternate 
impulsive dosage-response relationships. 


Cbmp*tiKl(iy with Single and Multiple Family 
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Cwtighted Day- 
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1.7% 

50 dB 

46 dB 

57.4 

Marginally compatible with tingle family, 
extensive outdoor use 

33 

55 

51 

58 

Marginally compatible with multiple family, 
moderate outdoor use and with multi -story, 
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63 

60 

55.6 

59.4 

Compatible with insulated multi-story use; 
incompatible with single and multiple family 
use 

123 

65 

60.2 

61.9 

Incompatible with any residential land use 

363 

75 

69.4 
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% Highly Annoyed - 100 / [1 + exp(lU3 • 0.141L in )] (FICON, 1992) 

% Highly Annoyed - 100 / [1 + erp( 11.17- (Galloway, 1981) 

% Highly Annoyed - 2.324 LCdn - 131.6 (Oklahoma City data) 
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are inherently more difficult to treat than those of airport neighborhoods. Sonic 
boom exposure produced by military operations is generally sporadic rather than 
regular, highly variable and difficult to predict accurately (due to the vagaries of 
long range acoustic propagation and uncertain flight tracks), and likely to be 
associated with nonlinear physical effects. Regular exposure to sonic booms near 
future overland flight corridors may be somewhat easier to predict. 

Not only are there more complications and loose ends in dealing with impulsive 
exposure than with general transportation noise, but there are far fewer data about 
exposure effects. Although several DoD and NASA-sponsored laboratory and 
field studies on the annoyance of sonic booms have begun to contribute new 
information, the body of information available for analysis is still considerably 
smaller than for the case of general transportation noise. 

Laboratory work, although useful for understanding how individuals’ immediate 
annoyance is affected by various aspects of impulsive signals, does not directly 
produce the sorts of information needed to generate and interpret impacts at the 
community level. Controlled field studies of longer term individual reactions can 
serve as a bridge between laboratory and community studies, but are difficult to 
design and conduct in ways that can test basic assumptions about the applicability 
of the equal energy hypothesis to relatively infrequent sonic booms. 

New technology, new assumptions, and new analyses are needed to identify and 
test improved means of predicting the effects of sonic booms on exposed 
individuals and communities. 


Table III. Relationship between numbers of sonic booms and land use 
compatibility. (Calculations performed for N-waves with 0.5 ms rise time and a 
duration of 350 ms, as described by Shepherd and Sullivan, 1991.) 
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ENDNOTES 


1. The term exposure is commonly used in two ways. One use of the term implies 
the time integral of intensity, while the other use implies the average sound intensity 
over a specific time period. Intensity is the rate of flow of sound energy per unit area 
per second. At distances from sound sources that are of interest in environmental 
analyses, sound intensity is directly proportional to the square of sound pressure. 
Thus, sound exposure is usually represented as the time integral of squared sound 
pressure. This process is often referred to informally as "energy summation". 
Magnitudes are reported in logarithmic terms. For example, sound exposure level is 
10 times the logarithm to the base 10 of the ratio of sound exposure to a reference 
exposure of 400 pPa -seconds. In this logarithmic form, squared sound pressure is 
called sound level and expressed in units of decibels. Sound level in decibel notation 
is often expressed as an average (equivalent) sound level over a specified time interval 
(usually 1 hour or 24 hours). Single events are often described by their sound 
exposure level (SEL) with a reference time interval of one second. 

2. This definition of "land use compatibility" does not deal directly with effects of 
noise exposure on people. Furthermore, certain of the remedies commonly used in 
airport neighborhoods to treat incompatible land uses (e.g., re-zoning land, insulating 
residences, purchasing avigation easements) are inappropriate in the case of en route 
exposure to sonic booms. Fidell (1992) contains additional discussion of the standard 
approach to dealing with the compatibility of aircraft noise and land use. 

3. According to Galloway (1981), "High-energy impulsive sounds of concern for 
community response are ... those for which the C-weighted sound exposure level... in 
any 2-second period is greater than 85 decibels (or greater than 75 decibels at night) 
and is 10 decibels greater than the C-weighted sound exposure level due to other 
sources in any contiguous 2-second period. These levels correspond to peak 
overpressures greater than approximately 105 decibels (95 decibels at night), that is, 
greater than approximately 0.1 pounds per square foot. 

4. The C-weighting network was selected in lieu of the A-weighting network in a 
noise metric intended to characterize high energy impulsive noise because two 
impulsive sounds with very different low frequency energy content may have the same 
A-weighted sound pressure level. This follows equally from the insensitivity of the A- 
weighting network to energy at frequencies below about 50 Hz, and from the fact that 

the spectral peaks of common impulsive noises are often two octaves yet lower in 
frequency. 

Approximate equivalences between C-weighted measurements and A-weighted 
measurements of a class of sounds (such as sonic booms) may nonetheless be 
established in several ways. One such equivalence (between A-weighted and C- 
weighted SEL values for sonic booms) is as follows: 


CSEL = 0.68(SEL) + 40.5 dB 

This relationship is derived from A- and C-weighted measurements of actual sonic 
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booms made outdoors (Galloway, 1981) and of recordings of simulated sonic 
booms made indoors (Pearsons, Tabachnick, Howe, Ahuja, and Stevens, 1993). It 
is appropriate primarily for estimating A-weighted levels for characteristic sonic 
booms with well-formed N shapes, but may also work reasonably well for sonic 
booms which have propagated long distances through the atmosphere, or for sonic 
booms occurring near the lateral cutoff distance. 

5. It was tacitly accepted in preserving Schultz’s assumptions that for purposes of 
predicting community response due to impulsive noise exposure, startle was fully 
accounted for by annoyance. It was likewise tacitly assumed that annoyance 
associated with secondary emissions was fully accounted for by substitution of the C- 
weighting network for the A-weighting network. ("Secondary emissions" are indoor 
rattling sounds produced by nonlinear re-radiation of low frequency impulsive energy 
from household contents.) 
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Atmospheric Turbulence 

The propagation of sonic booms through the real atmosphere can have a pronounce cITcct on 
t lie signature received on the ground. It has been well established that turbulence in the losscr 
part of the planetary boundary layer known as the mixing layer is a significant contributor to 
the distortion of sonic booms, as illustrated in this figure. The changes in the atmospheric 
conditions during a day and from day to day results in a large variation in the sonic boom 
signature measured on the ground for an aircraft (lying at a nominal operating condition at 
different times oT the day. I he objective of this study is to evaluate the variability in the 
loudness of the booms due to these propagation effects. 


ATMOSPHERIC TURBULENCE EFFECT 



192 


BOOM FILE Database Descriptions 


The BOOM FI LI! database contains overpressure distributions for 43 passes of the several 
types of aircraft previously mentioned. This data was collected using Boom Even >_ 

Recorders (BEAR) and modified Larson Davis LD700 Personal Dosimeters. I hose dev ecs were 
arranged in a linear array of 13 microphones located perpendicular to the mglu path at sideline 
distances ranging from 0 miles (i.c., directly under the flight path) to roughly 2 mi . ■ ■ 

in the figure. g Thc aircraft fiew across the microphone array with steady (light conditions which 
were achieved several miles prior to reaching the microphones. BOOM! ILL also contains 
aircraft tracking data which consists of altitude, Mach number, climb angle, acceleration, 
heading, and lateral and longitudinal position with respect to a reference microphone I his lata 
is provided at one second intervals Tor most of the aircraft overnights. Limited atmospheric 
data was also collected during the BOOM FILL tests. This data consisted of ground station 
wind speed and direction, air pressure, and air temperature measured just prior to each set o 
flvovcrs. Upper atmosphere rawinsondc data recorded at nearby weather stations on ic . 
days are also provided. This consists or wind speed and direction, sound speed, relative 
humidity, dew point, temperature and pressure at 1,000 foot altitude intervals ranging ron 
roughly 2,500 to 100,000 feet above mean sea level (Reference I). 


BOOMFILE MICROPHONE ARRAY 
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BOOM FILE Flight Conditions 

The BOOM FILE flight conditions arc listed in this table. The range of conditions is large from 
Mach 1.0 to 3.0 and altitude 13,000 ft to 70,000 ft with 9 different aircraft. However, repeat runs 
of the same aircraft at similar flight conditions arc limited. 


BOOMFILE Flight Conditions Summary 


FLIGHT TRACK MACH ALTITUDE BOOH AT SITE 00 
DATE AIRCRAFT INTERSECTION HUMBER (Ft HSLT (Local Time) 


31 JUL 87 

F-4 * 

57.8 

1.20 

16000 

08:41:20 

03 AUG 87 

F-4 

60.1 

1.24 

29200 

07:48:33 


F-4 

60.6 

1.29 

29300 

07:58:33 


F-4 

53.6 

1.10 

13000 

08:08:04 


F-4 

59.2 

1.10 

14400 

10:29:59 


F-4 

61.3 

1.37 

44400 

10:43:22 


T-38 

58.6 

1.00 

13600 

10:05:35 


T-38 

56.0 

1.10 

13000 

10:12:15 


T-38 

59.5 

1.11 

29600 

12:28:18 


T-38 

60.5 

1.05 

21200 

12:38:17 

04 AUG 87 

AT -38 

60.0 

.1.17 

41400 

07:19:41 


AT -38 

60.0 

1.12 

32300 

07:30:09 


AT -38 

63.0 

1.15 

16700 

07:36:46 


AT-38 

59.6 

1.20 

30300 

09:14:06 


AT -38 

59.0 

1.10 

14000 

09:23:15 


F- 15 

61.5 

1.38 

41400 

07:56:42 


F- 15 

60.3 

1.20 

29700 

08:04:06 


F- 15 

60.6 

1.10 

12500 

08:10:13 


F- 15 

60.0 

1.13 

15200 

10:46:15 


F- 15 

59.0 

1.28 

31000 

11:02:18 


F- 15 

64.0 

1.42 

45000 

11:11:28 


F- 15 

60.0 

1.40 

45500 

11:34:21 

05 AUG 87 

F- 16 

57.0 

1.25 

29500 

09:06:05 


F- 16 

60.0 

1.43 

46700 

09:33:54 


F- 16 

58.8 

1.17 

19300 

09:44:51 


F- 16 

59.5 

1.13 

14400 

11:44:24 


F- 16 

60.6 

1.12 

13800 

11:54:39 


F- 16 

60.5 

1.25 

30000 

12:04:46 


SR-71 

60.8 

2.50 

64800 

09:26:12 


SR - 71 * 

59.8 

3.00 

73000 

10:55:12 


SR * 71 

59.4 

1.23 

32400 

11:08:38 


SR - 71 

62.0 

1.70 

52000 

12:35:51 

06 AUG 87 

F- 18 

60.0 

1.30 

30000 

07:44:12 


F- 18 

59.6 

1.40 

44700 

07:57:05 


F- 18 

58.0 

1.10 

14200 

08:10:36 


F- 18 

59.8 

1.30 

30000 

10:22:47 


F* 18 

59.8 

1.43 

45000 

10:34:14 


F- 18 * 

59.8 

1.10 

13000 

10:48:38 


F- 14 

56.2 

1.20 

31500 

08:28:45 

r 

F- 14 

62.0 

1.27 

16500 

10:43:43 


F- 1110 

59.8 

1.20 

14000 

11:48:18 


F-111D 

59.8 

1.40 

45000 

12:04:44 

07 AUG 87 

F- 1110 

58.3 

1.25 

29900 

10:50:26 


For each of these flights, except where noted by an 
asterisk, tracking data are provided 
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XB-70 Database Descriptions 


The XB-70 database (Reference 2) consists of frequency spectra and overpressure time histories 
Tor 39 flights (51 runs) of the XB-70 aircraft. The data was collected at several ground stations 
using a microphone, tuning unit, d.c. amplifier, and I'M tape recorder setup played back into 
a recording oscillograph. I he oscillograph plots were then digitized using an optical scanning 
system. In this test program the microphones were arranged in one of two configurations, 
either three ground and one pole or six ground and two pole microphones all located within a 
200 foot by 200 foot grid pattern shown in the figure. The location of the measurement site 
with respect to the aircraft flight path for different runs ranged from directly overhead to a 
sideline distance of over 15 miles. Each run is considered as one flight over one cluster of 4 or 
8 microphones. Atmospheric data for the XB-70 database consists of digitized trace plots for 
temperature and wind speed parallel and perpendicular to the flight path for all runs. I he 
National Oceanic and Atmospheric Association provided pressure, temperature, wind, and 
relative humidity vs altitude profile rawinsondc data at 12:00 and 24:00 hours. They also 
provided limited test site climatological data consisting of temperature, wind speed and 
direction, cloud cover description, and dew point within an hour of each run. This database 
has more repeat runs than BOOM FI IXi, however, the sideline distance to the microphone 
cluster varied significantly from run to run (Reference 2). 


XB-70 MICROPHONE CLUSTER ARRANGEMENTS 
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XB-70 Flight Conditions 


The flight conditions of the XU-70 database arc listed in this table. 


XB-70 SONIC ROOM LOG 

(for flights of March 4, 1965 through May 27, 1966) 


DJM 


5/Cl- 

T/O 

T/O 

Fit. 

Boos 

Boon 

Boon 

Boon 

Land 

riui 

Date 

Fit 1 

Time 

Gr.Wt 

. Tlse 

Tise 

Mach 

Alt 

Gr.Mt. 

(Jr. Wt 

i 

3-4-68 

1-7 

1018 

480K 

1:37 

iiin 

1.83 

50500 

3371 

297 X 

2 

4-20-65 

1-10 

1113 

510K 

1:42 

1213 

1.00 

48000 

3501 

300 X 

3 

7— 1-65 

1-14 

0650 

610K 

1:44 

0600 

2.60 

66000 

310K 

205K 

4 

7—2 7—o s 

1-15 

0707 

510K 

1:43 

0732 

1.23 

32000 

4231 

300X 

5 

8-10-65 

2-2 

0700 

470X 

1:27 

0740 

1.38 

42300 

357K 

31QX 

S 

8-18-65 

2-3 

1220 

490 X 

1:58 

1330 

1.40 

46000 

38 IK 

305X 

1 

8-20-65 

2-4 

1115 

493K 

2:04 

1159 

1.42 

42500 

38 7 K 

2B5K 

8 

9-22-65 

1-16 

1200 

510X 

1:57 

1225 

1.50 

33800 

456JC 

300X 

9 

9-29-65 

2-6 

1147 

495K 

2:04 

1220 

1.35 

33000 

4 40 1C 

295K 

10 

10-5-65 

2-7 

1213 

495C 

1:40 

1243 

1.42 

31000 

438K 

2B5X 

11 

10— 11-65 

2-8 

1310 

515K 

1:55 

1332 

1.51 

34000 

4231 

298 X 

12 

10-14-65 

1-17 

0906 

510K 

1:47 

0938 

1.76 

41000 

433K 

300X 

13 

10— 16-65 

2-9 

0912 

520K 

1:43 

1027 

1.40 

50000 

313K 

295K 

14 

11-2-65 

2-11 

1126 

520K 

1:54 

1255 

1.80 

50500 

317K 

295X 

15 

11-4-65 

1-18 

1019 

5151 

2:04 

1105 

1.87 

41500 

357 X 

300 X 

l b 

11-18-65 

1-21 

1233 

515E 

2:02 

1338 

1.61 

41500 

348X 

300X 

n 

1 1-30-65 

1-22 

0000 

515K 

1:59 

1010 

1.02 

53000 

325X 

295X 

18 

12-1-65 

2=13 

0902 

525K 

2:02 

1030 

2.31 

60000 

328X 

297K 

1 9 

12-2-65 

1-23 

0925 

516K 

1:50 

1040 

1.79 

54000 

317X 

300X 

20 

12-3-65 

2-14 

0906 

5201 

1:55 

1030 

2.48 

65500 

329K 

300K 

21 

12-10-65 

1-25 

1230 

515S 

2:18 

1315 

1.55 

30500 

436K 





(2nd 

run) 

— — 

1400 

1.25 

38000 

371X 

295X 

22 

12-11-65 

2-15 

oasa 

520K 

2:03 

0918 

1.50 

37000 

454K 





(2nd 

run) 

— 

1020 

2.90 

70000 

321K 

300 X 

23 

12-21-65 

2-16 

1307 

510K 

1:49 

1427 

2.92 

70000 

32 IK 

300X 

24 

1-3-66 

2-17 

0901 

520K 

1:52 

1020 

2.91 

69800 

317K 

295K 

25 

1-11-66 

1-31 

0702 

447K 

1:35 

0750 

1.80 

44900 

369K 

2B5K * 

26 

1-12-66 

2-18 

0855 

525K 

1:46 

1018 

2.05 

86000 

297 X 

290E 

27 

1-15-66 

1-33 

HOB 

450K 

1:27 

1153 

1.78 

45100 

373X 

290 K 

26 

3-4-80 

1-38 

1055 

523K 

2:27 

1140 

1.75 

41000 

446X 



(2nd 

station-ease run) 



1140 

1.82 

42000 

445X 

293X 

29 

3—7—66 

1-37 

1402 

5201 

2:19 

1532 

1.17 

41000 

344X 



(2nd 

etatlon-easo run) 


1532 

1.17 

40000 

343K 

295K 

30 

3-15-66 

2-24 

0809 

5351 

1:59 

1030 

2.66 

68500 

310X 



(2nd 

etatlon-eose run) 

— 

1030 

2.66 

C9300 

310X 

293X 

31 

3—17—66 

2-25 

0847 

5351 

1:52 

1015 

2.74 

66000 

308 K 



(2nd 

etatlon-ees4 

\ run) 

— 

1015 

2.74 

66000 

308 K 

297X 

32 

3-10-66 

2-26 

1040 

530 JC 

1:57 

1210 

2.84 

70300 

30 5 X 



(2nd 

etatlon-eani 

i run) 



1210 

2.04 

70300 

304 X 

291X 

33 

3-28-60 

1-40 

0050 

520K 

1:41 

1053 

1.00 

51000 

319K 



(2nd 

atatlon-ease 

» run ) 

— - 

1053 

1.00 

51000 

319K 

300X 

34 

3-29-66 

2-29 

1027 

530 JC 

1:51 

1137 

1.56 

44000 

314K 



(2nd 

etat ion-ease 

i run) 

— — 

1137 

1.56 

44000 

314K 





— (2nd 

run) 

— 

1152 

1.36 

36400 

30 4 X 



(2nd 

eta tion-2nd 

run ) 

— 

1152 

1.36 

36400 

304 K 

300 X 

35 

4-5-60 

1-42 

1026 

520K 

2:01 

1138 

1.55 

52000 

334K 

295K 

36 

4-21-66 

1-45 

1539 

524K 

2:02 

1646 

2.26 

63000 

338X 

290X 

37 

4-23-66 

2-35 

1120 

525JC 

2:01 

1140 

1.11 

32000 

468K 



(2nd 

etetlon-ease 

run) 


1140 

la 18 

32000 

467K 

_ _ 





— (2nd 

run) 

— * 

1255 

2.20 

64000 

362X 



(2nd 

etation-2nd 

run! 

— — 

1255 

2.20 

84000 

362K 

310X 

38 

5-16-66 

2-38 

0900 

520K 

2:00 

1040 

1.30 

44300 

321K 

300 X 

39 

5—27—66 

2-42 

1100 

520K 

2:08 

1240 

1.24 

39000 

310K 

300K 


Total number of sonic boom flights = 39 


Total number of sonic boom ; 

runs 

r-4 

fl 
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Analysis Progression 

This figure shows the organization for the remainder of the presentation, beginning with a 
description of the variables calculated in what will be referred to as the extended database. The 
analysis of these calculated parameters with respect to the aircraft flight conditions and flight 
times is reviewed. The analysis then is focused on time of day variations. This is followed by 
further analysis in terms of lateral cutoff and morning vs afternoon comparisons. 


OUTLINE 


• EXTENDED DATABASE 

• VARIABILITY IN OVERPRESSURE, RISE TIME, AND LOUDNESS WITH 
FLIGHT CONDITIONS 

• SONIC BOOM VARIABILITY IN REPEAT FLIGHTS 

• VARIABILITY WITH TIME OF DAY 

• VARIABILITY IN BOOM SYMMETRY 

• STATISTICAL DISTRIBUTION BY LATERAL DISTANCE 

• STATISTICAL DISTRIBUTION BY TIME OF DAY 
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Extended Database 


Several noise metrics and various categories of rise time were determined for each point in the 
BOOM I'll J! and XB-70 databases. These quantities, listed in the figure, arc available in tabular 
and digital format. An example of this extended database is shown. 


ANALYSIS PARAMETERS 

• MAXIMUM OVERPRESSURE, P max 

• STEVENS MARK VII PERCEIVED NOISE LEVEL, PLdB 

• A-WEIGHTED SOUND EXPOSURE LEVEL, ASEL 

• C-WEIGHTED SOUND EXPOSURE LEVEL, CSEL 

• RISE TIME FROM 10% TO 90% P max 

• RISE TIME TO 50% P max 

• RISE TIME TO 75% P max 

• RISE TIME TO 100% P max 
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Prediction Method 


I he noise metrics for the measured boom signatures were compared to metrics for signatures 
predicted using Carlson’s simplified method (Reference 3) option of the sonic boom analysis 
program MDBOOM (Reference 4). In this technique, a simple F-function input is sealed to the 
local conditions. I he scaling factors used are the lilt parameter, Ki . determined from the aircraft 
Mach number, weight, length, and local pressure, and the shape parameter, K s determined 
from the aircraft type and K, as shown in the figure. K s is then used to scale the simple 
F‘ -function of the figure by the factor shown, flic signature is then evolved to the microphone 
(far field), resulting in a change of amplitude. An aging or steepening calculation is then 
performed to arrive at the signature propagated through a non-turbulcnt atmosphere (ideal 
N-wavc). 


SIMPLIFIED SONIC BOOM PREDICTION PROCEDURE 



Enter lift parameter 
Vm 2 - 1 W 
Kl " 1.4p v M 2 t 2 
Select shape factor Kg 
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Flight Condition Groups 

For sonic boom variability analysis the BOOM FI I-H and XB70 data were each divided into 
four groups based on aircraft altitude and Mach number values. The range of flight conditions 
for these groups arc shown in the figure. 


FLIGHT CONDITIONS 

BOOMFILE DATABASE 


RANGE 

GROUP 1 

GROUP 2 

GROUP 3 

GROUP 4 

Altitude (feet) 

~ ■ 

10,000 - 20,000 

25.000 - 35,000 

40.000 - 50,000 

50.100 - 80,000 

Mach number 

1.05 - 1.30 



o 

o 

s — 

1.10 - 1.50 

1.50 - 3.50 

Sideline Distance 
(feet) 

0 - 45.000 

0 - 55,000 

0 - 80.000 

0 - 60,000 


XB70 DATABASE 


RANGE 

GROUP 1 

GROUP 2 

GROUP 3 

GROUP 4 

Altitude (feet) 

30,000 - 40,000 

40.100 - 50,000 

50.100 - 60.000 

60.100 - 72,000 

Mach number 

1.17 * 1.55 

1.17 - 1.87 

1.55 - 2.31 

2.05 - 2.92 

Sideline Distance 
(feet) 

0 - 50,000 

0 - 80,000 

0 - 70,000 

0 - 80.000 
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Overpressure Variability Dependence on Flight Conditions 

Comparisons between measured and calculated data were made Tor the various metrics and rise 
times as they varied with sideline distance. I lie high altitude / high Mach number group, 
bottom figure, shows very good agreement between measured maximum overpressures and 
predicted maximum overpressures (assuming uniform noil-turbulent atmosphere). For the low 
altitude / Mach number group, upper figure, the comparison is relatively poor. While the 
measurements in both groups include the cfTccts of propagation through the lower layer of 
turbulent atmosphere, the high altitude / high Macli number group represents flights where 
measurements arc well within the lateral cutoff, where the boom has already propagated an 
adequate distance so that the 'shock is in an equilibrium state prior to entering the lowest .5,000 
Tcct of the atmosphere, and where the propagation through the atmosphere is more vertical than 
the low altitude / low Mach number flights. These can be expected to reduce variability in 
measurements and improve theory - data agreement. These plots include morning as well as 
afternoon flights. 

OVERPRESSURE VARIABILITY 



0 5 10 15 20 25 30 35 40 45 


SIDELINE DISTANCE (X1000 FT1 




□ = MEASURED 
• = PREDICTED 


AIRCRAFT TYPE ALL 
MACH NUMBER RANGE 1.70 TO 3.00 
ALTITUDE RANGE (FT) 50000 TO 80000 
FLIGHT TIME 6:00 TO 14:00 
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Rise Time Variability Dependence on Flight Conditions 

The variability in the rise times (measured from 10% to 90% P m „) for two groups of 
measurements is plotted in these figures. Again, the low altitude / low Mach number group (top 
figure) shows a wider range of values (up to .'>0.3 msec) compared to the smaller variation (up 
to 1 1.8 sec) for the high altitude / high Mach number group (bottom figure). 
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Loudness Variability Dependence on Flight Conditions 

These figures illustrate that, as would be expected based on the observations of previous figures, 
the loudness of the measured and predicted booms arc in good agreement for the high altitude 
/ high Mach number flight groups (bottom figure). I or the low altitude / low Mach number 
altitude group the loudness values of the measured booms arc scattered around the predicted 
boom loudness values (top figure). For the other frequency domain metrics, ASliL and CSHL, 
similar comparisons were noted. 


LOUDNESS VARIABILITY 
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Sonic Boom Variability in Repeat Flights 

The BOOMl'ILH database lias a limited number of repeat flights (i.c., same aircraft at 
nominally the same flight conditions). Analysis of measurements from two such sets are shown. 
The first set consists of two flights of F16 aircraft at an altitude of 14,000 feet. I he second set 
consists of two flights of FI5 aircraft at a higher altitude of 45,000 feet, in each ease, the repeat 
flights were made within a few minutes of the first flight. These plots thus provide a real 
representation of the variability in sonic boom measurements due to propagation effects. It can 
be noticed that the lower altitude runs (top figure) show a greater variability in sonic boom 
maximum overpressure compared to the higher altitude runs (bottom figure). 


SONIC BOOM VARIABILITY IN REPEAT FLIGHTS 
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Sonic Boom Variability in Repeat Flights 


These plots show that the lower altitude runs (top figure) have greater variability rise time 
compared to the higher altitude runs (bottom figure). 
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Sonic Boom Variability in Repeat Flights 

These plm* show that the lower altitude runs (top figure) have greater variability in loudness 
compared to the higher altitude runs (bottom figure). 


SONIC BOOM VARIABILITY IN REPEAT FLIGHTS 



MACH ALTITUDE F15 
O = 1.42 45000. meas 
■ = 1.42 45000. pred 

* = i.40 45500. meas 

• = 1.40 45500. pred 



SIDELINE DISTANCE (X1000 FT) 
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Sonic Boom Variability in a Repeat Data Point 


Data from two different flights at similar flight conditions roughly 10 minutes apart show the 
effect propagation through the atmosphere has on sonic booms. The overpressure plots show 
that the 1 1:44 flight resulted maximum overpressure of 5. 19 psf and a rise time of 5 ms whereas 
the 1 1:54 flight yielded a maximum overpressure of 2.66 psr with a rise time of 6.875 ms. The 
difference in loudness between these two booms was .1.1 PldB. Ihc spectra plot shows that the 
11:44 flight had less low frequency noise (below 200 ID) and more high frequency noise (above 
200 1 1 7.) than the 1 1:54 flight. 


SONIC BOOM VARIABILITY IN A REPEAT DATA POINT 


Mach 1.13 Altitude 1M00 ft Sideline 9.500 ft 
Mach 1.12 Altitude 13.800 ft Sideline 9.600 ft 



Frequency, Hz 
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Variability with Time of Day 


The measurement'; set up in the XB-70 database used either a three or six ground microphone 
cluster in a 200 by 200 foot square. Only minor variations arc expected from one microphone 
to the other in the absence of significant propagation effects. Atmospheric turbulence and thus 
the propagation arc expected to vary with the time of the day. The figure examines the 
variation in maximum overpressure with time of day. The data points arc for flight conditions 
Mach = 1.17 to 1.87 and altitude = 40,000 (4 to 50,000 0 (identified as Group 2 previously). 
The variation in values front one cluster to another is due to differences in operating conditions 
and sideline distances but the variations within a cluster are due to propagation differences. It 
can be noticed in this figure that the variability in maximum overpressure is very small for 
morning flights (prior to 1 1 AM). Around noon and in the afternoon this variability increases 
a little. 
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Variability with Time of Day 

This figure shows the variation in rise time with time of day. I he data points are for flight 
conditions Mach = 1.17 to 1.87 and altitude = 40,000 ft to 50,000 0 (identified as Group 2 
previously) The variation in values from one cluster to another is due to differences in 
operating conditions and sideline distances hut the variations within a cluster arc due to 
propagation differences. The rise time shows a large increase in variability in the afternoon. 
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Variability with Time of Day 


This figure examines the variation in loudness with time of day. The data points arc for flight 
conditions Mach — 1.17 to 1.87 and altitude — 40,000 ft to 50,000 ft (identified as Group 2 
previously). The variation in values from one cluster to another is due to differences in 
operating conditions and sideline distances hut the variations within a cluster arc due to 
propagation differences. It can be noticed in this figure that the variability in loudness can be 
as much as 10 PI.dB in the measured data. Similar variability in loudness was noticed in groups 
1,5, and 4 of the XB-70 database with the higher altitude runs generally having slightly lower 
variability. 
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Variability in Boom Symmetry 

Sonic boom asymmetry was determined by the difference between overpressure, PLdB, ASH I,, 
or ('SUL calculated separately Tor the compression portion and the expansion portion of the 
sonic boom signature. The variability in A overpressure Tor the lower altitude / lower Mach 
number group of flights (top figure) is slightly greater than the high altitude / high Mach 
number group or nights (bottom figure). The lower values and smaller variability in A 
overpressure (front shock) for the higher altitude / higher Mach number group is consistent with 
the near N-wavc signatures and reduced atmospheric effects associated with these signatures in 
this altitude group. 
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Variability in Boom Symmetry 


Sonic boom asymmetry was determined by the difference between overpressure, PLdB, ASEE, 
or ( SEE calculated separately for the compression portion and the expansion portion of the 
sonic boom signature. The variability in A PEdB in the lower altitude / lower Mach number 
group of flights (top figure) is slightly greater than the high altitude / high Mach number group 
of flights (bottom figure). In the "afternoon hours", the asymmetry in loudness has a greater 
variability than the asymmetry in overpressure (previously shown). This is an indication of the 
larger effect of atmospheric turbulence on rise time. 
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Statistical Distribution 


The XB-70 database was divided into two data groups - tliosc within .50 percent or the 
calculated lateral cutofT distance (dye) and those outside of this boundary. Such a grouping has 
been used in Reference .5 in the analysis of BOOM FI F.H data. The histograms these figures 
represent the distribution of measured maximum overpressure values (normalized by the 
corresponding calculated uniform atmospheric maximum overpressure) in the database for these 
two groups. It can be seen that for the below 50% dye group (bottom figure) maximum 
overpressure distribution is nearly symmetric and has approximately a normal distribution 
shape with shorter skirts (smaller variance) away from the mean value. By comparison, the 
above 50% dye group (top figure) shows a large variability in measured maximum overpressure. 
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Statistical Distribution 


'I lie XB-70 database was divided into two data groups - those within 50 percent of the 
calculated lateral cutoff distance (dye) and those outside of this boundary. Such a grouping has 
been used in Reference 5 in the analysis ofBOOMIILn data. The histograms in these figures 
represent the distribution of the difference between measured and calculated (uniform 
atmosphere assumed) loudness in the XB-70 database Tor these two groups. The below 50% 
dye group (bottom figure) has a symmetric distribution with a -0.15 dB mean for PL„„ S - PL„ k 
whereas the above 50% dye group (top figure) has a bi-modal type distribution with a -1.7 dB 
mean and larger variance about the mean. 
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Statistical Distribution by Time of Day 

The variability of measured maximum overpressure in the below 50% dye group was further 
analyzed in terms of the time of day to statistically quantify the turbulence effects. The 
histogram in top figure shows that the maximum overpressure measurements for the morning 
(before Ham) flights have a smaller variability than for flights which occur aflcr 1 lam (bottom 
figure). While the mean values of maximum overpressure in the two plots arc not very different, 
the mean values occurs more frequently before I lam than aflcr 11am (i.c., the afternoon 
distribution has a larger variance). This trend was also observed in the sonic boom 
measurement program at White Sands Missile Range (Reference 6). 

STATISTICAL DISTRIBUTION BY TIME OF DAY 
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Statistical Distribution by Time of Day 

The variability of measured loudness in the below 50% dye group was further analyzed in terms 
of the time of day to statistically quantify the turbulence effects. The histogram in top figure 
shows that the loudness measurements for the morning (before I lam) flights have a smaller 
variability than for flights which occur after I lam (bottom figure). 
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Summary 


The BOOM FILE and XB70 databases were analyzed in various calculated metrics and rise 
times which are available in electronic format. The variation in boom loudness was observed 
to be as much as 10 PTdB with larger variation occurring in the lower altitude / lower Mach 
number flights, afternoon flights, and outside the carpet semi-span. Analysis of asymmetry 
showed that differences of up to 12 dB occurred and was greater for lower altitude / lower Mach 
number flights. 


SUMMARY 

• ANALYZED VARIABILITY IN SONIC BOOM SIGNATURE PARAMETERS AND 
METRICS (DUE TO ATMOSPHERIC PROPAGATION) USING BOOMFILE AND 
XB-70 DATABASES 

• 10 dB VARIABILITY IN BOOM LOUDNESS POSSIBLE 

• GREATER VARIABILITY OBSERVED IN 

LOWER ALTITUDE / LOWER MACH NUMBER FUGHTS 
AFTER MID-MORNING FLIGHTS 

MEASURMENTS OUTSIDE 50% OF SONIC BOOM CARPET SEMI-SPAN 

• UP TO 12 dB ASYMMETRY OBSERVED 

• VARIABILITY IN BOOM ASYMMETRY GREATER FOR LOWER ALTITUDE / 
LOWER MACH NUMBER FLIGHTS 
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Figure 1 

Outline of presentation. 


SONIC BOOMS AND MARINE MAMMALS 

Bill Cummings 
Oceanographic Consultants 
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Figure 2 


Examples of baleen (blue whale, top) and toothed (bottlenose dolphin, 
bottom) whales. 


Balaenoptera musculus 
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Figure 5 

Estimated population sizes of selected large whales (from Gaskin). 


WHALE 

POPULATION SIZE 

Sperm Whale 

over 567,800 

Bowhoad Whalo 

about 2,500 

Right Whale 

about 3,500 

Gray Whalo 

between 7,000-15,000 

Humpback Whale 

between 3,000-5,000 

Blue Whale 

between 7,000-13,000 

Fin Whale 

about 68,000 

Sel Whalo 

about 130,000 

Bryde's Whale 

over 60,000 ■ 

Mlnke Whale 

between 1 13,000-646,780 


The best estimates of large whale population sizes (data from 
Gaskim 1982). 


224 



Figure 6 


Sound source levels of selected species of marine mammals (various 
sources). 


MARINE MAMMAL SOUND SOURCE LEVELS 

dB re 1 ^iPa. 1 m 

Baleen Whales 

131 - 189 

Toothed Whales 

135-200 

Seals 

95-180 

Sea Lion, Walrus 

157-172 
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SPL (DECIBELS re 1 uPa) 


Figure 7 

Audiograms of five toothed whale species (from indicated sources). 


AUDIOGRAMS OF FIVE TOOTHED WHALES 
(adapted from Fay, 1988) 



1- Bottlenose dolphin (Johnson) 

2- Amazon River dolphin (Jacobs & Hall) 

3- Killer whale (Hall & Johnson) 

4- Harbor porpoise (Andersen) 

5- Beluga whale (White et al.) 
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SPL (DECIBELS re 1 jiPa) 


Figure 8 

Audiograms of five pinnipeds (from indicated sources). 


AUDIOGRAMS OF FIVE PINNIPEDS 
(adapted from Fay, 1988) 



1- California sea lion (Schusterman et al M 1972) 

2- Harbor seal (Mohl, 1968) 

3- Ringed seal (Terhune & Ronald, 1975) 

4- Harp seal (Terhune & Ronald, 1972) 

5- Northern fur seal (Schusterman & Moore, 1978) 
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Figure 9 


Sound spectra from a sonic boom at the water’s surface and at depths 
indicated, compared to spectra of typical, low and high oceanic ambient 
noise (after Cook). 


PENETRATION OF A SONIC BOOM INTO THE OCEAN 
{adapted from Cook et al., 1972) 



FREQUENCY (Hz) 
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SPL (dB RE 1 jiPa) 


Figure 10 


Sound spectrum of a sonic boom as measured at 15 ft depth, compared to 
average hearing curves of selected pinniped and toothed-whale species. 


Sonic Boom Spectrum Vs Pinniped 
and Tooth ed-Whale Hearing Curves 
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Figure 11 


This study showed that marine mammal distributions will coincide with 
SS aircraft sonic booms. Depending upon the species and location, 
marine mammals will sense, not sense, or questionably sense sonic 
booms. Marine mammals probably will not be harmed physically. It is 
not probable that whale or swimming pinniped behavior will be 
significantly affected by sonic booms. Hauled out pinnipeds may or may 
not be affected by booms possibly depending upon species, age, and 
behavior at the time of the boom. There is little doubt that marine 
mammals will habituate to some extent. 


COINCIDENCE + 

SENSATION + - ? 

DIRECT PHYSICAL HARM - 

BEHAVIORAL EFFECTS 
Whales - 

Swimming Pinnipeds - 
Hauled-Out Pinnipeds + - ? 

HABITUATION + ? 
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Figure 12 

There has only been one experimental study (Sea World) on one species 
of pinniped using simulated sonic booms and none based upon real 
booms. No studies have been undertaken on the effects of sonic booms 
on any other marine mammals. I recommend early contacts with the 
National Research Council and Acoustical Society of America 
committees, with the regulatory authorities at the National Marine 
Fisheries Service, and with the Marine Mammal Commission. Low cost, 
comparatively easily facilitated studies may be undertaken on captive 
porpoises and pinnipeds using simulated booms for the purpose of 
determining any effects or habituation. However, regulatory 
authorities may not require such studies. A marine mammal permit 
should be applied for well before any decision is made regarding 
technical feasibility of the overall program. 


EARLY CONTACTS (NRC, ASA, NMFS, MM COMMISSION 

CONDUCT SYMPOSIUM (CONJUNCTION WITH SMM) 

SSB TESTS ON CAPTIVE PORPOISES & PINNIPEDS 
(ONLY IF REQUIRED BY NMFS) 

HABITUATION TESTS 
(ONLY IF REQUIRED BY NMFS) 
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